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ABSTRACT 



This dissertation is an integrated method of designing sampled- 
data feedback control systems to required time response. This method 
is distinguished from the classical frec^uency design methods in two 
ways - first, the closed-loop rather than the open-loop function serves 
as the starting point of the method; and, second, the procedure invol- 
ves nearly direct synthesis rather than trial and error. To implement 
the design procedures, a set of correlation theorems are derived which 
permit interpretation of the closed-loop system time responses in 
terms of the pole-zero configuration of the closed-loop system function. 
Further, a graphical method is presented which permits rapid solu- 
tion of the open-loop poles when the closed-loop poles are known. The 
techniques and results of an analog computer study are presented which 
serve as an experimental verification of the methods developed in this 
dissertation. 
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INTRODUCTION 



!♦ 1 General Background 

The past few years, since the end of World War II, has witnessed 
a tremendous increase in industrial and military applications of feed- 
back control systems. As a result of this activity, the field of feed- 
back control systems has developed and specialized in many direc- 
tions. One important class of systems that has recently been sub- 
jected to vigorous investigation is known as sampled-data systems ( or 
sampling servomechanisms). In this class of systems one or more .of 
the variables consists of trains of pulses or sequences of values. These 
pulses are generally available at equally spaced intervals of time and 
no information is received between two consecutive pulses. Conven- 
tional systems, or systems in which the variables are always continu- 
ous functions of time, may, in fact, be considered as a special limit- 
ing case of sampled-data control systems, where the sampling interval, 
or period, approaches zero. However, the present body of published 
sampled-data control system theory is most useful when applied to sys- 
tems that have sampling rates low compared to the maximum response 
frequency of the system. It is to that type of system that the findings 
presented in this manuscript are directly applicable. 

One of the earliest applications of sampled-data systems was a 
tracking system including a search radar. In this case the antenna, 
rotating at a constant angular velocity, results in information about 
the location of the object which is available only as the antenna sweeps 
past the object. Time multiplexing of telemetering data is another com- 
mon source of sampled-data. The increasing trend toward systems 
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containing a digital computer as one element is furnishing further 
examples. 

In the aforementioned applications the processes require that the 
information be in sample form. In other cases sampling is purposely 
introduced although the signal information is available as a continuous 
time variable. Generally this is done for the express purpose of 
increasing sensitivity or for more economical use of equipment. The 
electromechanical galvanometer and chopper bar is a classical exam- 
ple of such a device (1), 

Among the early contributors to the study of these systems were 
Oldenbourg and Sartorius (2), They applied difference equations to the 
analysis of sampled-data control systems. The solutions of these equa- 
tions yielded the absolute stability of the system as well as the response 
to various inputs. This approach offered the same difficulties as exper- 
ienced in the classical methods of using linear differential eqxiations 
for describing continuous signal feedback control systems. The method 
is cumbersome and does not permit relating the performance of the 
systems to its constants. As a result of these difficulties, attention 
has been given to developing transform calculus methods similar to 
those available for continuous systems. 

Through the independent efforts of Hurewicz (3), Stibitz and Shan- 
non (4), transform methods were evolved which gave input and output 

Ts 

relations of the system as a function of z, where z is defined as e 
(s is the Laplace complex variable. ) This actually is a special case 
of the Laplace transform in which the transform is a rational function 
of z . Because of this, the inverse transform yields the time solu- 
tion only at the sampling instants. 
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Ragazzine and Zadeh (1), pursuing the developments of Hurewicz, 
formalized a transform calculus now known as the z-transform. Trans- 
form pairs similar to Laplace transform pairs were nnade available. 
Transfer functions in terms of the variable z were shown to be deriva- 
ble in straight forward manner. Even more important, they demon- 
strated how sampled-data could be analyzed by z-transform methods. 
This article laid a firm foundation for future frequency analysis inves- 
tigation. Significant among the many recent contributors are Jury (5) 

(6) (7) (8) and Lago (9) (10). They showed that methods and concepts 
associated with the analysis of continuous systems could with some 
modification be directly applicable to sampled-data control systems 
using the z-transform. Correlation was established between the fre- 
quency response and the transient response of some second order sys- 
tems. This infornnation was shown to be useable for approximate design 
of higher order systems where dominance holds. 

Linvill (11) showed that the transfer function in a sampled-data sys- 
tem could be represented in the frequency domain by an infinite series. 
Fortunately, only two or three terms are generally needed for a good 
approximation. The labor involved is about the same as with z- trans- 
form techniques but it has the advantage of being easily applied by 
those familiar with the standard methods of analyzing continuous vari- 
able systems. 

As the previous discussion indicated, most investigators of sam- 
pled-data control systems have concerned themselves with analysis. 
They have shown that either the z-transform or Laplace transform 
methods can be applied to the open- loop system and readily yield infor- 
mation as to the response of the closed- loop system. Although the 
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labor involved is considerable greater than for continuous systems of 
the same complexity, both methods are straight forward. That the 
methods are applied to the open-loop system is significant, since in 
general this is known in advance by the designer. Furthermore, com- 
pensating networks, employed to improve the closed-loop responses, 
are added as part of the open-loop system. To date, design procedures 
are largely based on analysis and trial and error methods. Predic- 
tion of high order system time domain performance is based on study 
of second order systems. In these respects the design of sampled- 
data systems is similar to the commonly employed methods for design- 
ing continuous closed- loop systems. 

In continuous system design considerable use is made of the abi- 
lity to subdivide the system into transfer function blocks, the require- 
ment'being that each block not affect the operation of the adjacent blocks. 
This requirement is satisfied by many physical devices and consequently 
the open-loop functions of continuous servo systems can readily be rep- 
resented by tandem transfer functions. An analogous situation exists 
with respect to sampled- data control systems. However, in this case 
the transfer function must be applied to all components between output 
and input of samplers. If a system employs only one sampler then one 
z-transform function must be used to describe the entire open-loop 

A 

system. Whenever a new network, such as compensating network, not 
employing a sampler, is introduced into the system an entirely new 
transform must be recalculated for the complete system. Even the 
slightest change in a parameter value of an existing component will 
require the determination of a new z-transform. (Fortunately, the 
introduction of a sampler into an otherwise linear system does not 
affect linearity; consequently, gain enters into the z-transform as a 
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multiplying factor. ) Because of these complications, the locus shap- 
ing on the Nyquist diagram or Bode diagram requires prodigious amounts 
of labor when compared with the amount of work required for a simi- 
lar analysis of continuous systems. If this complication could be removed 
design of sampled-data control systems would be no more difficult than 
for continuous ones. 

This investigator spent considerable effort studying this problem 
in hopes of finding some design methods which would circumvent or 
ameliorate this complication. It readily became apparent that the solu- 
tion was not to be found by changes or additions to present methods 
which modify the open- loop transfer function. Further progress could 
only be had by adopting a completely new approach. 

Since working with the open-loop transfer function necessitates 
considerable labor and complications, it was thought that perhaps the 
closed-loop transfer function might offer opportunities for design with- 
out the proportional increase in labor experienced when this approach 
is used in continuous systems (12). Since system performance is gen- 
erally specified in terms of closed-loop performance, this approach 
would actually have the advantage of being a design method based on 
systhesis as contrasted with the present methods based on analysis. 
Investigation along this line produced a direct synthesis method which 
permits the designer to shape the time response of the closed-loop 
system. Also, the procedure indicates and helps to resolve possible 
conflict between the performance specifications and the plant (fixed 
components in the system). When applied to higher order systems, 
this method requires less time than those methods whose starting point 
is the plant transfer function. The detailed results of this investigation 

\ 
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constitutes the contents of this dissertation. 



1,2 Discussion of the Dissertation 

After a few introductory comments the subject material of this dis- 
sertation will be discussed by chapters. The figures showing system 
time responses were determined by an analog computer. These com- 
puter curves were found to check out the analytical calculations for 
values at the sampling instants exceedingly well, the error, in general, 
being less then reading errors of the 4 centimeter paper. So far as 
possible, the nomenclature used in this dissertation is that recom- 
mended by the AIEE Standard Subcommittee on Terminology and Nomen- 
clature of the Feedback Control Committee (13). Because of the nature 
of this dissertation additional terminology was required. 

Chapter II presents various aspects of current sampled-data theory 
which will be required for subsequent developments in this disserta- 
tion. Because of the newness of the field, important developments are 
found scattered throughout the literature and it would not be expected 
that the average reader have available all the references from which 
the material of this chapter is drawn. Hence, Chapter II presents this 
material in more detail than would be normally expected. For the rea- 
der familiar with the field, this chapter could be scanned or omitted 
completely. No new contribution is made in this chapter but the author's 
explanations may add further insight into the theory. 

Chapter III is the major contribution of this dissertation, A philo- 
sophy of sampled-data system design is developed which is nearly pure 
systhesis. With this method the designer is able to take time domain 
specifications and proceed in a logical manner to the final system. 
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The method permits including specific components as part of the spe- 
cifications. 

The procedure is quasi graphical leaning strongly on root locus 
techniques (14), This will find favor with practicing engineers who 
generally prefer graphical over analytical methods since they are less 
subject to mistakes and usually provide more circumspection. 

The procedure may be described by the following five steps: 

(1) The poles and zeros of the closed-loop function are deter- 
mined from the time domain specifications, 

(2) The poles of the open- loop function are found by graphical 
methods, 

(3) The open-loop poles and zeros are shifted to include the 
important poles and zeros of the specified components. 

(4) The new closed- loop poles and zeros are determined gra- 
phically and checked against the specifications, 

(5) The appropriate compensation function is determined. 

Steps (3) and (4) may be eliminated if the designer feels that 

the system will operate satisfactorily with a compensator cancelling 
out all the poles and zeros of the specified components. 

Chapter IV develops the necessary techniques to implement the 
design philosophy described in Chapter IIL Step (1) requires a set of 
correlation theorems which permits interpretation of the closed-loop 
system time responses, such as peak output, steady state error, set- 
tling time, etc, in terms of the pole-zero configuration of the closed- 
loop system frequency function. Some correlation theorems are avail- 
able in the literature (6) (15), However, the set of error coefficient 
theorems and the settling instant theorem presented in this chapter 
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were derived by the author. 

The terms closed-loop root locus and open- loop root locus have 
been introduced to distinguish between the former which is the well 
known "Evans Technique" (14) and the latter which is a method of 
finding the open- loop poles from the closed- loop function. In this dis- 
sertation the open- loop root locus was a necessary development to 
implement step (2). This technique, originated by the author, repre- 
sents an important addition to the present root locus methods for sam- 
pled-data control systems. Step (3) presents no difficulties while step 
(4) utilizes the standard root locus techniques. Current literature has 
numerous instances where the closed-loop root locus or Evans locus 
has been applied to sampled-data control systems (6) (12) (15), Only 
portions of these references which are directly applicable to this design 
method are included. The details of step (5) are left to Chapter V, 

Basic pole-zero configurations are derived which give the engin- 
eer a starting point for his design. It is shown that most specifica- 
tions can be satisfied by a basic pole-zero configuration consisting of 
one pair of complex dominant poles plus several additional sets of real 
poles and zeros which must be carefully located according to the rules 
developed in this dissertation. 

Chapter V applies the technique developed in Chapters III and 
IV to a wide variety of sampled-data control systems which included 
hold networks and delays. Although not an exhaustive list, the success 
of the design methods on these problems indicates its general applica- 
bility, The steps involved in deriving the Laplace transform from the 
z-transform are illustrated. To the author's knowledge this is the 
first time that this process has been presented in such detail. The 
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experimental findings verifying the methods of this dissertation are 
presented. 

The general theory is extended to types of sampled-data control 
systems which do not enjoy a closed- loop transfer function independent 
of the input function. This is a further complication found in sampled- 
data systems not having a direct analog in continuous linear systems. 
These systems are physically characterized by the omission of a sam- 
pler in the error channel. Conventional analysis methods readily test 
the absolute stability. However, to the author's knowledge this disser- 
tation is the first to include design procedure for these very special 
types of control systems. 

Methods for handling systems with multiple inputs, such as a 
disturbance input, are also discussed. 

Chapter VI presents the analog computer set-up that was used 
to verify the methods in this thesis. The use of an analog computer, 
with the transfer function of the various components making up the sys- 
tem directly analoged, permits the system to be analyzed without the 
effects of nonlinearity overshadowing the features of the theory to be 
tested in the experiment. A specially designed mechanical sampler 
( 16 ) and several new experimental techniques are discussed. 
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INTRODUCTION TO SAMPLED-DATA THEORY 
Z, I Operation of the Sampler 

A linear sampled-data feedback system can be represented by the 
introduction of a switch in the otherwise continuous linear system, as 
shown in Figure 2. 1. This switch, referred to as the sampler, is the 
component that distinguishes a sampled-data system from a continu- 
ous system. Figure 2. I is misleading regarding the operation of the 
sampler, for in physical systems the sampler operates in accordance 
with the circuit shown in Figure 2. 2 . Here the switch, SW, is alter- 
nately closed on position a for T seconds and closed on position b 
for seconds. The travel time between positions a and b is con- 
sidered negligible and the sampling period T becomes 

T = (2.1) 

In this dissertation the sampling period is assumed constant. If more 
than one sampler is present in the system, they are considered as 
operating with the same period and synchronized. 

With the sampler illustrated in Figure 2,2 , the output signal is 
other than zero only for the seconds each period the switch is in 
position a , Thus the operation of the sampler converts a continuous 
signal e(t) at its input into an output signal of pulses with width T 

a 

and a height proportional to the amplitude of e(t) at the sampling instant. 
Figure 2. 3 illustrates the input and output signals of a physical sam- 
pler. 

In the systems to be considered in this dissertation, the time T 

a 

is much smaller than both T^ and the time constants of the systems. 
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Figure 2. 1. A sampled-data system 





Figure 2. 2. Physical circuit for sampler 
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(a) Continuous signal applied to the input of the 
sampler 
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(b) Pulse train at output of the sampler 
Figure 2. 3. Input and output signals of a physical sampler 
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Physically this latter requirement means that the system responds a 
negligibly small amount while the pulse is applied. If this is the case, 
the system response would be pie same to an impulse with strength 
equal to the area of the pulse. Since mathematically the impulse is 
much simpler to handle than the pulse, both in the literature and in 
this dissertation, the impulse train is substituted for the actual value 
of the continuous time variable. It should be noted, that the area of 
the impulse is considered as equal to. the value of the continuous vari- 
able at the instant of sampling rather than corresponding to the area 
of the pulse. With the area of the pulse equal to the value of the con- 
tinuous variable multiplied by T , the mathematical representation 

3l 

of the sampler involves a gain change. Thus, the gain K of the open- 
loop system is related to the actual gain used in the physical sys- 
tem by: 




The Figure 2. 4 may yield more insight to the nature of the approx- 
imation involved when the pulse is replaced by an impulse under the con- 
ditions stated above. It is clearly shown in this figure that the fre- 
quency content of the pulse is nearly constant over the extent of the 
bandwidth of the open- loop function. Thus, the system output contains 
essentially the same proportion of frequency components with either 
the pulse or impulse input. The change in gain is illustrated with T 

sl 

less than one second. 

It is helpful for analysis to consider the sampler as an impulse 
generator which is modxilated by the input signal. The process of con- 
sidering sampling as modulation of an impulse train is shown in Fig- 
ure 2.5 , wherein the values of the impulses are represented by the 
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(a) The pulse and impulse time functions 




(b) The frequency spectrums of the impulse (1) , pulse 
(2) and bandwidth of the open- loop system (3) 



Figure 2, 4. The impulse approximation illustrated in terms 
of the Frequency Domain, 
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(b) Impulse train 




Figure 2, 5 . The sampler as an impulse modvilator 
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the amplitude of the lines occurring at the sampling instants. With 
input and output denoted as e(t) and e*(t) respectively, the output 
may be written as 



e*(t) = e(t) i(t) 

where i(t) represents a train of unit impulses 



(2.3) 



i(t) 



I “oO-nT) 



(2.4) 



m 



Equation 2. 3 may be written equivalently as 



e*(t) = ^ e(nT) u^ (t - nT) 



(2.5) 



Here the negative values of n are omitted by virtue of the assump- 
tion that e(t) is zero for all negative values of time. The Laplace 
transformation of both sides of Equation 2. 5 yields 



An alternate and very useful expression for E*(s) can be 
obtained by the use of contour integration. Equation 2.3 represents 
the sampled time function as the product of two time functions from 
which the Laplace transform can be obtained by a process of complex 
convolution. The complex convolution integral is given by (17) 



oo 




( 2 . 6 ) 



n=o 



L[£, WfjO] = ^ 



/ 



c + 

Fj(s - w) (w) dw 



(2.7) 



c - j oo 
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The contour of the integration of Equation 2. 7 is a path from c - j ©o 
to c + j'*’ along an analytic strip between the singularities of 
Fj (s - w) and those of (w) . 

Letting the impulse train be represented by f^ (t) of Equation 
2. 7 , Fj (s - w) becomes 

0-0 

F(s - w) a I(s - w) s ^ e" (2.8) 

nao 

5 

This infinite summation can be expressed in closed form as 

I (s - w) a (2,9) 

1 - 

Then, equating e(t) to f^ (t) yields 

F^ (w) a E (w) (2.10) 



With the indicated substitutions in Equation 2.7 the Laplace transform 
of e*(t) becomes 



L [e* (t) ] a 



1 



2irj 




1 



dw 



1 - e 



T(s-w) 



( 2 . 11 ) 



C - J OO 

Evaluation of Equation 2. 11 can be effected by closing the poles 
I (s-w) in the right half of the w-plane as shown inUgure 2, 6 , When 
the degree of the denominator of E(s) in s is higher than that of the 
numerator by at least two, the integral is zero along the infinite semi- 
circle. The alternate path in the left hand plane will be considered in 
a later section of this chapter. 

Applying Cauchy's integral formula gives 
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Figure 2.6, Paths of integration used to evaluate the 
complex-convolution integral. 
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z?3“^ 



E(w) 



1 



— dw = Res 

1-e poles of l_ 

I{s- w) 



E(w ) 
n 



, -T(s-w ) 
1-e n J 



( 2 . 12 ) 



The poles of the I(s - w) are the zeros of 



1 . e- T{s-w) ^ 0 



(2. 13) 



or 



l j . . 2 TT n 

u/ = s + J 

n T 



(2.14) 



These poles, though infinite in number, are all simple and hence the 

residues are given by Expression 2. 15 
n = oo 



E(w) 



z - 

d/dw [1 . 



(2.15) 



Since 



d/dw 






T(s- w) 



^ . . 2irn 

oyss-t-j -Tj .— 



Te 



j 2 tr n 



( 2 . 16 ) 



The value of the integral represented by Equation 2. 12 is 



oo 



E«(s) ^ yY. ^ 



, 2 IT n 



) 



n=- oo 



(2. 17) 



Defining the sampling frequency as LJ where 



H = z,/ 



(2.18) 



Equation 2, 17 may be written as 



19 






M 



n= <yo 



E*(s) 



T n=- (X) 
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^ E(s + j n(Jg) 



(2. 19)* 



It is seen that two equivalent expressions are available for analy- 
tical description of the sampler, as demonstrated by Equation 2. 16 
and Equation 2, 19, The first equation, which expresses the Laplace 
transform of the sampler output in terms of the values of the input 
function at the sampling instants, is particularly helpful in gaining 
insight into the physical operation of the sampler and its effect on sys- 
tem performance. The second equation places in evidence the fre- 
quency effect of the sampler on the input signal. 

Replacing s by s + j , where m is any integer, in Equa- 

tion 2, 19 or its alternate expression Equation 2, 6, results in an 
expression which is identical with E* (s) (1) (18) , That is, E* (s) 

is a periodic function of s with period , Expressed symbo- 

s 

lically 



This indicates that the s-plane can be divided into strips of equal width 

with centers at n(^ as shown in Figure 2, 7 , 

s 

The spectrum alteration performed by the sampler is shown in 
Figure 2. 8. From this it is seen that the sampler contains an exact 
replica of the input as well as an infinite number of additional images 

♦Equation 2, 19 is in error if e(t) 0 at t = 0+ , In this case the 
Equation 2, 19 should be modified to read 



E* (s + j miJJ - E* (s) 

s 



( 2 . 20 ) 



n= oo 




(s + j n U ) + 1/2 e(0+) 

s 



n= - oo 
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(a) Amplitude spectrum of E(s) 




(b) Amplitude spectrum of E* (s) 



Figure 2. 8. Input and output amplitude spectrum of 
a sampler. 
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centered about a multiple of the sampling frequency. These additional 
components are referred to as complementary signals (11) . Spaces 
between images would indicate that it is possible to regain the exact 
signal by filtering. Unfortunately, a satisfactory filter for the process 
would introduce too much time delay. Any sigfial with frequencies 
greater than the sampling frequency would overlap into the adjacent 
strip with the result that the information would be lost. This conclu- 
sion is equivalent to Shannon's sampling theorem which requires that 
the sampling frequency be at least twice the highest input frequency. 



2,2 Elements of z-Transform Theory 

If the path of integration for Equation 2. 11 had been closed to 
the left, the finite number of poles of E(w) would have been enclosed 
and there would be a finite number of residues to evaluate. This alter- 
nate path is also indicated in Figure 2.6, In this case 



L [e* (t) ] s ^ Res 

poles of 
E(w) 

The use of Eqxiation 2.21 may be illustrated by determining the 
z -transform of the sampled sequence of a unit ramp, 

e(t) . = t (2.22) 



E(w ) 
n 



-T(s-w ) 
e ' n' 



( 2 . 21 ) 



The addition is a result of the difference in definition of the sam- 
pling process in Equation 2. 6 which implies a full sampling at t'*0 
and Equation 2. 11 which assumes the impules to be an even function 
sampling only one- half of the value e(o) . In feedback control sys- 
tems Equation 2, 19 is normally adequate. 
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The Laplace transform is then 

E(s) = (2.23) 

s 

E(s) is seen to have a second order pole at origin. Performing the 
operations indicated by Equation 2.21 yields 






E* (s) = 



T e 
sT 



sT 



(2.24) 



(e“" - 1) 

An alternate method of deriving Equation 2.24 employs Equation 2.6 
as a starting point. Hence 



e*(t) s Tu^(t - T) + 2Tu^(t - 2T) + 3Tu^(t - 3T) + . . . (2. 25) 

The Laplace transform of this series of samples is 



E(s) s Te‘^® + 2Te"^^® + 3Te"^^® + ... (2.26) 



Multiplication of both sides of the above equation by yields 



E* (s) ^ I ^ __L_ + J_ ^ _ 

„ Ts 2Ts 3Ts 4Ts 

T e e e e 



(2.27) 



sT 

Multiplication of both sides of the above equation by de and 
integration gives 



f de®^ = 

J Te®^ 



I 

2Ts 

e 



-i ... + C (2.28) 

3Ts 



where C is the constant of integration. The series on the right of 
the above equation maybe rewritten in closed form. 
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E* (s) , sT 

i- ' d e = 

T 



sT , 
e - 1 



+ C 



(2. 29) 



sT 

Differentiation of Eqxiation 2.29 with respect to e yields 



T e®^ 

E* (s) = (2.30 

- l)^ 

sT 

which is identical to Equation 2.24. The substitution z = e into 
the expressions gives rise to the z-transform. The above expression 
in z-transform notation becomes 



E(z) 



T z 




(2.31) 



3 

By these methods it is possible to derive a table of z-transforms 
(5). Table 2. 1 is a short list of z-transforms. 

The inversion of a z-transform function to obtain the time res- 
ponse of the function at sampling instants can be effected by any^one 
of three ways: 

(1) Evaluation of the inversion integral (19) 



c(n t) 




(2.32) 



where P is a path of integration enclosing all the singular 
points of C(z) 

(2) Partial fraction expansion of the fiinction so that each term 
can be identified with the corresponding time function of 
Table 2. 1. 

(3) Series expansion in descending powers of z in which the 
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Table 2. 1. Short List of z- Transforms 
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coefficient of corresponds to the value of the tinne 

function at the nth sampling instant. 



2.3 z-Transfer Function ' ' 

Of primary importance in servo analysis are the input-output 
relations of networks. This relation for samipled data systems will 
now be established using the simple systems illustrated in Figure 2.9i 
which shows a linear network W pulsed by a sampler. In addition 
to the desired continuous output time signal a fictitious sampler is pro- 
vided to generate a sampled output signal. Although not actually part 
of the signal flow in the system, consideration of the sampled output, 
rather than the continuous output, simplified the analysis as will be 
shown. 

It is evident from Figure 2. 9 that the Laplace transform is 
given by 



where G(s) is the Laplace transform of network W and E* (s) 
is the sampled input. Applying Eqxiation 2.19 yields the transform 
of the pulse output of the output sampler. 



Recalling the periodicity of pulsed-transfer function as expressed by 
Equation 2.20 , the above expression may be rewritten 



C(z) s G(s) E* (s) 



(2.33 




G(s + j n (J ) E* (s + j n U) (2. 34) 

o S 



C* (s) 




E* (s) G(s + j n ) 

8 



(2.35) 
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Figure 2.9. A pulsed network 
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Figure 2. 10. A system with a sampler between 
networks. 
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Since E* (s) is a factor in each term of the above series, it may be 
placed before the summation, hence 



C* (s) s E* (s) I G(s + j n L/g 

n= -oo 

Denoting the bracketed term by G* (s) 

G* (s) = - 1 ^ ^ G(s + j n 

ns- oc 



(2. 36) 



(2. 37) 



Equation 2. 36 may be rewritten 

C* (s) s E* (s) G* (s) (2.38) 

In the z- transform notation 

C(z) = E(z)G(z) (2.39) 

It is seen that the z-transform of the output of a network is simply the 
z-transform of the input multiplied by the z-transfer function. Fur- 
thermore, the z-transfer function is derived in the same manner as 
the transform of a time function. 

An alternate expression for G*(s) in terms of g(t) , the 
impulse response of W , can be derived by noting the similarity in 
form between Equation 2. 19 and Equation 2.37 given below for com- 
parison 

oo 

E* (s) = -\r- Y. E(s + j n U^) 

ns- CX5 

oo 

G* (S) = Y J ” ^g) 

ns- oo 
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Taking the Laplace transform of Equation 2. 3 yields 



E* (s) = L [e(t) \ i(t) ] (2.40) 

Since G{s) is the Laplace transform of g(t) it follows at once that 

G* (s) = L [ g(t) X i(t) ] (2.41) 

The over-all transfer function of two networks separated by sam- 
plers is readily seen to be equal to the product of the z-transform for 
each network. This is illustrated in Figure 2. 10 . If no sampling 
device separates the networks, the input to the second network is the 
continuous function of time present at the output of the first network. 
Then the corresponding pulsed (starred) transfer function will be 
I n= 

G* (s) = -4--^ G^s + j n (v'^) G^(s + j n CJJ (2.42) 

na-oc 

Hence the over-all transfer function must be calculated by considering 
the system as a whole. The accepted notation (1) for this case is 

G* (s) a Gj G *2 (s) (2.43) 

or in terms of z-transform notation 

G (z) a Gj G^ (z) (2.44) 

The transfer function of a system with feedback, such as shown 
in Figure 2. 11 may be derived in the following manner. The expres- 
sion for the error is 

E(s) a R(s) -H(8 )C(s) (2.45) 

The output is directly related to the sampled error by 
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Figure 2. 11. A simple feedback sampled-data system 



I 



31 



C(s) S E* (s) G(s) 



(2.46) 



Equation 2,46 substituted into Equation 2.45 yields 

E(s) = R(s) -H(s)G(s)E* (s) 

Taking the starred transfer function of both sides 

E* (s) = R* (s) -HG* (s)E* (s) 



Solving for E* (s) 



E* (s) 



The output is then 



C* (s) 



R* (s) 

1 + HG* (s) 

R* (s) G(s) 
1 + HG* (s) 



Taking the starred transfer function of both sides 



C* (s) 



R* (s) G* (s) 
1 + HG* (s) 



and in z-transform notation 



R(z) G(z) 
1 + HG(z) 



(2.47) 



(2.48) 



(2.49) 



(2. 50) 



(2.51) 



(2. 52) 



Table 2.2 gives the transforms of the output for a number of 
basic servo configurations. Table entries 2, 6 and 7 show the input 
may not always appear as a simple multiplicative factor in the expres- 
sion for C(z). In this case it is not always possible to derive an 
expression for the feedback system transfer function independent of 
the input. As discussed in the introduction this creates a type of sys- 
tem whose analog is not found in linear systems operating wholely on 
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Table 2. 2. Output Transforms for Basic Sampled-Data Systems 
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continuous data. 



2,4 Final and Initial Value Theorems, Error Coefficients and Mapping 
The output of a system C(z) may be expanded in partial fractions 
as 

D z ^ 1 ^ 

C(z) = + Jpj. + (2,53) 

z - 1 1 2 

z-e z-e 

The inverse transform into the time domain is readily determined by 
entry 4 in Table 2. 1 

-d,T -d,T 

c(t) s + D^e ^ + D^e +... (2.54) 

Assuming that the real parts of the d's are positive all terms of Equa- 
tion 2. 54 go to zero as time goes to infinity with the exception of . 

Therefore, the final value of c(t) may be found directly from C(z) 
as follows 

lim(z - 1) C(z) = Urn c(t) (2.55) 

z -^1 t 

Expanding C(z) in descending powers of z gives 

H H 

C(z) = H + + — ^ + ... (2.56) 

O 

z z 

The inverse transform into the time domain is 

c(t) s H +H,U(t-T) + H,U^(t - 2T) + ... (2.57) 

O 1 O C 0| 

Hence, , the initial value of c(t) , maybe found directly from 

C(z) by 
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(2.58) 






lim C(z) - lim c(t) 
z — ^oo t -^0 

The error expression for single loop systems with unity feed- 
back is 

R(z) 



E(z) = 



For a step input 



Thus 



R(z) = 



E(z) = 



1 + G(z) 



Z - 1 



(z-1) [1 + G(z) ] 



(2.59) 



(2.60) 



(2.61) 



Applying the final value theorem the steady state error e is 

s s 



, . z-1 

e = Um 

ss z 

t — &0 



(z-1) (1 + G(z) ) 



= lim 



z —^1 1 + G(z) 



(2.62) 



Then the error coefficient for a sampled-data system may be defined 
as 



K* 



or 



Controlled Variable 
Actuating Error 

1 - e 



K* 



s s 



s s 



= lim G(z) 
z 



In like manner it can be shown 



(2.63) 



(2. 64) 



K* 



lim (z - 1 ) G(z) 
z -^1 



and 



K* 



1 

T 



lim (z-1) G(z) 
z — 



(2.65) 



( 2 . 66 ) 



It was explained in Section 2. 1 that the s -plane could be consid- 
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ered as divided into strips. The conformal mapping of the first strip 
in the s-plane is the unit circle in the z- plane. This and mapping of 
other important lines are shown in Figure 2.12. 

2. 5 Response between the Sampling Instants 

The output of a sampled- data system is continuous even though 
its input is sampled, but the z-transform only contains information of 
the time response at the sampling instants. This section presents an 
extension of the z-transform called the modified z-transform which 
permits determining the output of a network between sampling instants. 

Consider the system in Figure 2. 13 which introduces a ficti- 
tious delay between the output of the network and the sampler at the 
output. A represents the fraction part of the delay to the total sam- 
pling period and hence has a positive value between zero and one. The 
effect of this delay on the signal presented to the sampler is illustrated 
in Figure 2. 14 . For reasons which will become clear later it is aus- 
picious to define the change of variable 

A s 1 - m (2. 67) 

Equation 2. 67 implies that mT is a pure advance. A system incor- 
porating this pure advance before the sampler and a delay of one sam- 
pling period following the sampler, as shown in Figure 2. 15, is seen 
to have the same output as the system in Figure 2. 13, The impulse 
response of the network and pure advance is g(t + mT) where g(t) 
is the impulse response of the network alone. Employing Equation 
2,41 it follows: 

Cq (z) = L t g(t + mT) X i(t) ] (2.68) 
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the relationship z = e 
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Figure 2. 14. The output compared with the delay output of 
system in Figure 2. 11 . 
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Figure 2. 15. Sampled- data system with output equivalent to 
that of system in Figure 2. 13 . 
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Equation 2.68 may be evaluated using the complex convolution integral 
in the following form: 



Cq U) = 




C+j oo 



G(w) e 



mwT 



C-J CO 



1 - e 



\ 

- T(s-w) 



d w 



(2.69) 



Evaluating the integral in the left hand plane provides a closed 
form for (^) • should be noted that the substitution represented 

by Equation 2. 67 assured the convergence of the integral on the infin- 
ite semi-circle. 

Because output C^(t) is identical with (^) 6ut shifted in 

time by one period, it follows from the real translation theorem (17) 

(s) = e"®'^C*^(s) (2.70) 

Varying the value of m from zero to unity, the entire response 
between the sampling instants can be found. The z-transform evolved 
from this technique is called the modified z-transform. With this nota- 
tion the transfer function of the system in Figure 2. 13 , including the 
linear network and delay AT , is 



G (z, m) 




Residue 



poles of 
G(s) 



G(s) e 

, sT 
1 - e 



msT 



-1 

z 



(2.71) 



Table 2.3 is a short list of modified z-transforms for various 
forms of G(s) . 



2. 6 Analysis and Design of Sampled- Data Control System by 
Conventional Techniques 

The absolute stability of sampled-data control systems can be 
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Table 2, 3, Short List of Modified z-Transforma 
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determined in much the same manner as is done with continuous sys- 
tems, For continuous systems a very powerful technique has been 
devised involving the construction of asymptotic diagrams in which 
log modulus and phase are plotted against log frequency. This is 
known as the Bode diagram and is primarily useful when the transfer 
function is expressed as a product of simple factors and the variable 
increases indefinitely in magnitude. That latter requirement is not 
characteristic of the z-transform, where for this diagram, z is 
varied along the unit circle. 

A Nyquist diagram may be constructed by plotting the gain and 
phase of the open-loop transfer function as z is varied along the unit 
circle. Closure of the coordinate - 1 + jO indicated instability. Con- 
struction of appropriate circles (20) gives an indication of the 

relative stability. 

As an example of analysis and design using the Nyquist diagram, 
consider the system shown in Figure 2, 16, The sampling period is 
0. 5 seconds. The first step in the analysis involves determining the 
open- loop z-transfer function, G(z), This requires a partial fraction 
expansion of G(s) , converting each term to z-transform with the 
aid of Table 2. 1 and finally a collection of these terms over a common 
denominator. These steps are illustrated below. 



G(s) 


10 

S ■■ ■ 




(2. 72) 




s(s + 3.2) 




Q(s) 


3.125 


3. 125 


(2. 73) 


” 3 


s + 3. 2 


G(z) 


3.125 


3. 125 


(2. 74) 


z - 1 


z - ♦ ^ 
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O igure 2. 1 6. 



Samp led -data system to be compensated. 
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(2. 75) 



Analysis on the Nyquist plot. Figure 2. 17 shows the system to 
be unstable. Let the design require a peak overshoot of 50% in the 
time response while maintaining the same low frequency gain. As a 
first try, it will be assumed that this reqmrement could be met with 
phase lag network compensation 

s + 2 



^C2<«) = ♦ s .H - T73 - 



(2. 76) 



The new z-transform open-loop transfer function must be deter- 
mined as shown in Equation 2. 77 through 2. 80 . 



^ /.V _ 7. 5(s + 2) 

^1 ■ s(s + 1. 5)(s + 3.”2) 



Gj (8) = 



Gj (z) 



Gj (z) = 



3. 125 



3. 125z 
z - 1 



L.470 

s + 1T5 ■ 

1.47z 
z “57T 



1.655 
s + 3.2 

1. 655z 

z - . 2 



1. 78z(z - . 41) 
(z.l)(z-:472y(z'..2) 



(2. 77) 
(2. 78) 

(2. 79) 

(2.80) 



The resulting Nyquist diagram. Figure 2.17 , gives the locus 
touching the = 2. 5 circle. This indicates a time response over- 

shoot of 53% . Trying a second network compensator 



^c2 



= . 5 



s + 2 

TTT 



(2.81) 



The G(s) and G(z) of the new open- loop system are then 



G, (s) = 



.5(s -h 2) 
s(s + l )(s + 3. 2) 



(2. 82) 
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1. 57z(z - . 376) 

(z -■ HTz - :Z0W~.-7U 



(2.83) 






The resulting Nyquist plot. Figure 2. 17 , shows the locus to be 
outside the equal two circle, indicating a peak overshoot of less 

than 44% . Hence, the design is completed as far as this procedure 
is concerned. The first step in calculation of the time response at 
sampling instants is to form the output z-transforms as given in entry 
4 of Table 2. 2. This is done in two steps, first by forming the closed 
loop transfer function and then followed by multiplication of the trans- 
form of the input. The final result is 

C(z) = (2. 84) 

z'^ - 1.24z^ + ,58z^ - ,46z + .121 



The first few values of the time response at the sampling instants 
are placed in evidence by dividing the denominator into the numerator 
as 



C(z) 



1.57 

z 





(2.85) 



Although the actual peak overshoot is outside of specifications, the 
Nyquist procedure has given a system which at least has a reasonably 
good time response. 

The previous design problem has illustrated the definite trial and 
error nature of design of sampled-data control systems by conventional 
techniques. For higher order systems, the labor may be prohibitive 
for designing to fixed specifications. 

Root locus methods starting with open- loop transfer functions 
suffer from the same difficulties in determining the effect of a com- 
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pensation network. 



Examination of the open- loop transfer functions given by Equa- 
tions (2, 75) (2. 80) (2, 83) permits drawing four important conclusions 
about z-transfer functions derived from continuous components. 

(1) The 8-plane poles have direct images in the z-plane 

(2) The s-plane zeros do not have direct images in the z-plane 

(3) Additional z-plane zeros are introduced when taking the 
z -transform 

(4) In a z-transform the numerator is of one degree less than 
the denominator if the function has a zero value at zero 
time and a nonzero value at the first sampling instant. 
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Figure 3, 1, A compensated sampled-data control system 
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Figure 3, 2. A sampled-data control system employing a 
sampler with the compensator 
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Figure 3. 1. A compensated sampled-data control system 




Figure 3. 2. A sampled-data control system employing a 
sampler with the compensator 
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culties the classical approach has been to work directly with the open- 
loop function. A compensator is chosen to reshape the Bode diagram 
according to certain tenuous relationships between the frequency and 
time domain. For continuous control systems the uncertainty in the 
final performance of the closed-loop system is compensated for by the 
ease with which the method can be applied. 

The design example in Section 2.6 showed that this ease of pro- 
cedure for design based on the open- loop function does not apply to 
sampled systems unless a sampler is included between the compensa- 
tor and the plant. Such a case is shown in Figure 3. 2. It is these 
difficulties which lead to the philosophy of a design procedure with the 
closed- loop functions as a starting point. The utilization of this philo- 
sophy requires two new developments: 

(1) Method of determining the closed-loop transfer function 
from the specifications 

(2) Method of determining the corresponding open-loop trans- 
fer function. 

This dissertation evolves these two developments making possible a 
logical design procedure enabling the designer to take the specifica- 
tions and proceed in a direct path to the final system. 

3.2 Synthesis through the Closed- Loop Poles-Zeros Configuration 

The details of the synthesis procedure are shown by thfe flow 
chart of Figure 3.3. The various blocks represent steps in the pro- 
cedure. The letters beside each arrow identify the operation 
required to advance to the next block. 

The chart is headed with the time domain specifications. By 
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use of "correlation theorenns" it is possible to interpret these speci- 
fications in terms of closed- loop pole-zero configuration on the z-plane. 
From the pole-zero configuration the "open loop root locus" yields the 
corresponding open- loop pole-zero configuration. At this juncture the 
designer can choose either of two concluding procedures. 

(1) Proceed directly to determination of the compensator 

(2) Attempt to conciliate the differences between the poles 
and zeros of the open-loop function and those of the 
plant so that at least the nnajor poles and zeros of the 

( 

plant are included in the final open- loop function. 

If the first method is chosen the procedure is straight forward. 
The desired open-loop z-transform is converted into its equivalent 
Laplace transfer function. The transfer function of the compensator 
may now be found by simple division 



(s) 



Gp (s) 



(3.1) 



where (s) is the open-loop or "desired" transfer function cor- 

responding to the selected closed- loop function. The actual compen- 
sating network may be determined by standard network synthesis 
procedure (12) . 

A study of Equation 3. 1 reveals that the compensator will, in 
general, cancel out each pole and zero of the plant as well as add the 
poles and zeros of the desired function. In practice it is impossible 
to exactly cancel the poles and zeros of physical components such as 
electric and hydraulic motors since their precise characteristics are 
never known over the full range of operation. Fortunately, complete 
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Correlation theorems 
Open-loop root locus 
Inspect and shift 
Closed-loop root locus 
Convert and divide 
Network synthesis 



Figure 3.3. Synthesis Flowchart 
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cancellation is not necessary. System performance is not significantly 
changed if the pole and its cancelling zero only lie near each other. 

This situation gives ri^e to a pole and zero placed close together in the 
closed- loop function and the resulting residue of the transfer function 
in the pole is negligible. 

A more serious objection rests on the argument that it is not pos- 
sible to cancel out equipment time constants and expect the equipment 
to function in a linear mode over a large range of input values. It is 
unreasonable to expect that cancelling out a motor time constant and 
replacing it with a much smaller one will result in any increase in the 
actual maximum obtainable acceleration or that even the motor and 
gear combination will be able to follow a ramp which exceeds its max- 
imum speed. Hence for this reason it is strongly urged that the 
designer proceed according to method (Z) . 

In method (Z) the designer compares the poles and zeros of the 
desired open- loop with that of the plant. By shifting the location of 
the closed-loop poles-zeros it is usually possible to include the major 
poles and zeros of the plant into the desired open-loop transfer func- 
tion, The open and closed-loop root locus furnish a graphical method 
for bringing about this correspondence. The impos sibility of incor- 
porating the major components of the plant into the desired open-loop 
function does not imply failure of this design method. Actually, this 
is a signal to the designer that the fixed plant is not capable of pro- 
viding the necessary performance demanded by the specifications. 

This simply means that the plant must be improved or, if this is not 
practical, the specifications must be changed to concede the realities 
of the situation. Thus this method increases the designer's circum- 
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spection. 



If the compensator can be cascaded to the plant through a sam- 
pler as shown in Figure 3.2 , the preceding steps are still taken and 
the compensator can be found directly from the z-transform of the 
system. 



(2) 



^z) 
Gp (z) 



(3.2) 



3.3 Time Response and Basic Pole-Zero Configuration 

The critical step in this procedure involves choosing a pole- 
zero configuration to satisfy the time requirements. The typical 
time requirements for a transient response to a step is shown in Fig- 
ure 3.4. The steady state error to a ramp and parabolic input are 
shown in Figure 3. 5. Numerous pole-zero arrangements could be 
found to satisfy a particular set of time specifications. However, the 
correlations theorems to be formulated in the next chapter are based 
on a dominant p>air of complex poles which together with a real zero 
will fairly well control the transient response. A pole-zero spaced 
close together and near the point 1,0 in the z-plane can provide 
control over the steady state-error. With this, the basic pole-zero 
configuration in the z-plane is that shown in Figure 3. 6 . 

It is not to be inferred from Figure 3.6 that all specifications 
can be discharged by third order systems. However, in many cases 
the synthesis can be handled initially in this way. As a last step in 
the design additional poles or zeros may have to be added if the com- 
pensator is to be physically realized. (This is discussed in more 
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Figure 3.4. Typical transient response to a unit input 





(a) Ramp input 



(b) Parabolic input 



Figure 3. 5. 



Steady state errors to ramp or parabolic inputs 
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Figure 3 



F igure 3 
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detail in Chapter V. ) If the poles and zeros have large real parts in 
the s-plane they will give rise to z-plane poles and zeros near the ori- 
gin and have little effect on the previously determined results. Figure 
3. 7 shows a typical p>ole-zero configuration. 

It should be kept in mind that this control of time response 
extends only to the sampling instants. This generally is satisfactory. 
However, a complete description of the output between sampling instants 
may be secured by employing the modified z -transform technique to the 
final system. 

Plants with poles and zeros in the right hand half of the s-plane 
require special consideration. The poles will have images outside the 
unit circle in the z-plane and so will the zeros in most cases. Rather 
than attempt to cancel these poles and zeros they should be included 
as part of the final open- loop transfer function. 
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CHAPTER IV 



SYNTHESIS TECHNIQUES 



4, 1 Introduction 

This chapter develops techniques to implement the design phil- 
osophy described in Chapter III, First, simple relations between the 
pole-zero configuration and the time response are developed. The 
desired response, as shown in Figure 3, 4 and Figure 3, 5 , may be 
specified by the following five quantities: 

(1) The peak sampling instant, n^ , the sampling instant with 
the maximum overshoot for a step input, 

(2) The peak sampling overshoot, M , the amount of the 
first overshoot at the sampling instant n^ for a step input, 

(3) The settling sampling instant, n , the sampling instant 

s 

that the value of the output reaches within 5% of its final 
value for a step input. 

(4) The number of oscillation, N , for the interval up to the 

/ 

settling instant for a step input. 

(5) The steady- state error constant K* , K* , and K* 

^ p V a 

for a unit step, unit ramp and unit parabolic input respec- 
tively. 

The first four quantities describe the transient response while the last 
gives the steady- state error. 

Basic pole-zero configurations such as shown in Figure 3. 7 are 
next derived giving the engineer a starting point for his design. Fin- 
ally rules for applying the open- loop and the closed-loop root locus 
are presented. 
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4,2 Correlation Theorems 



For the sampled- data feedback control system shown in Figure 
4, 1 the closed-loop system function W(z) is 



W (z) = 



G(z) 

1 + HGrzT 



= K 



Q(z) 



(4.1) 



where 



Q(z) 

P(z) 



7T (z - qj^) 
k=l 



m 



Tf 

i=l 



(z 



Pi) 



(4.2) 



(4.3) 



The values of qj^ and p^ may be either real or complex. It 
is assumed that the expressions do not contain multiple roots and that 
the order of P(z) is greater than the order Q(z) . These two 
assumptions are true for most physical systems. When the input is 
a unit step, the output at the sampling instants becomes 



c(nTy s. Z'^ 




(4.4) 



where denotes the inverse transform. Thus employing 

Equation 2. 32 



c(nT) 




K Q(z) 

- 1) PU) 



z d z 



(4.5) 



where P is the contour of the integration in the z-plane that encloses 
all the singularities of the integrand in Equation 4, 5. Evaluation of 
the contour integral gives 
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Figure 4,2, Important characteristics of the output 
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where 



z 

o 



u + 
o 



j V 



o 



P'(s> 



dpfc) 

d z 



Z=Zq 



(4. ll) 



It is important to realize that the assumption in Equation 4.10 
does not neglect all other roots but neglects all other modes because 
these modes die out quickly. The accuracy of this assumption will be 
discussed later. 

In general, the peak of the output will not occur at a sampling 
instant. If n is permitted to be noninteger, then a value of n could 
be found, so that at one sampling period of time later the output will 
be repeated. This condition is illustrated by Figure 2.2 and can be 
expressed as 



c(n+l) - c(n) = 0 



(4. 12) 



The actual value of n^ is the upper integral value of the actual 
n obtained. In the derivation which follows it will be assumed that 
the solution of Equation 4. 12 results in an integer value for n . 

Substituting Equation 4.12 into 4.8 and solving for n yields^*' 



n 

P 




7T 

2 



ang Q(z^) + ang P(z^) 



(4. 13) 



^ The intermediate steps leading to Equation 4. 13 are quite involved. 
The more important steps are as follows: the constant term, the 
coefficient and the nth power of | | cancel yielding 
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Then the peak sampling instant can be expressed as 



0 n 
o p 



TT 



- (sum of angles from the zero to 
the predominate pole) - (sum of 
angles from the other poles to 
the predominate poles) (4. 14) 



cos 



n 0^ + ang Q(z^) - ang (z^ - 1) - ang P'(Zq) 



|z I cos 
‘ o ’ 



(n+l)Q^ + ang Q(z^) - ang (z^ - 1) - ang P*(z^) 



Expanding the right hand side of the above equation by the sum of 
angles trigonometric identity. 



cos n 0^ + ang Q(z^) - ang (z^ - 1) - ang P'(z^) 



IzqI 



n 0^ + ang Q(z^) - ang (z^ - 1) - ang f^(z^) cos 0^ 



- sin n 0^ + ang Q(z^) - ang (z^ - 1) - ang P*(z^) 



sin 0 



The following trigonometric relationships exist for the values above 



:os 0 



sin 0 



u 



|z I 
' o * 



cot 



ang (z^ - 1) 



u - 1 
o 



Substituting these into the preceding equation yields 
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Thus, it may be concluded that additional poles increase n and that 

P 

- 

a zero decreases n 

P 

The first overshoot is obtained by substituting Equation 4, 13 
into 4, 10 , Hence 



M = 



r/ 2.1/2 

[(u^-l) +v^ ] 



K Q(z ) 
P'<-o' 



|z I"p si 



sin 



(ang Zq- 1) I (4. 15) 



j (4. 



recognizing 



sin 



ang (z-1) 



[(u -1)^ + V 



1/2 



(4. 16) 



Eq\iation 4,15 simplifies to 



2 V 



M = 



[(u -1)^ + V ^ ] 

' o o ^ 



K Q(z^) 



p' 



1^0 l"p 



(4. 17) 



The accuracy of the assumption that all modes other than that 
due to the constant term and the predominate poles may be explained 
readily. From the complete solution, Equation 4.8, consider a cer- 
tain mode D(Pj^)^ . This mode will decrease to a small percent of 



Tan 



n 0^ + ang Q(z^) - ang P'(Zq) - ang (z^ - 1) 



Tan 



- ang (z^ - 1) 



From which Equation 4. 13 follows immediately. 
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its initial value when n = 3 , and the coefficient D is small if it 

P 

is not too close to the other roots. In the event that a certain mode 

must be accounted for the contribution of this mode at n can be 

P 

evaluated from the general Equation 4, 8, 

Unless the actual value of n from Equation 4, 13 is equal to an 
integer an error is introduced giving an optimistic value of M, This 
is because the value of the sine term in Equation 4, 15, for the non- 
integer n , is less than the sine term would be at n^ , This error 
is small unless the sampling period is a large proportion of the period 
of oscillation and n is near midway between two integer values. 

For a zero displacement error system the constant term must 

equal unity when the system responds to a unit step of position. Hence 

m 



n 

i=l (1 - p ) 

K = — L_ (4.18) 

k=l ^ 



Substituting Equation 4, 18 into 4. 17 the peak sampled overshoot can 
be expressed as 

M = (Product of distances from all poles to the point 
of 1,0), excluding distances from two predomin- 
ate poles to the point 1, 0) /(product of distances 
from all poles to predominate pole u + j v , 
excluding the distance between predominate 
poles) X (product of distances from zeros to the 
predominate pole u + j v )/(product of dis- 
tance from all zeros^ to th^ point I, 0) | z^ 

(4. 19) 

It is evident from Equation 4. 19 that the smaller n^ the 

larger M, Thus the choice of n^ and M requires a compromise. 

The effect of real poles on M is evident from Figure 4,3, For 

bj^ c 

position one near the origin — r; is nearly unity and the effect 

D a 
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Figure 4.4. 



The effect of a real zero on M and n 

P 
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of the pole is slight. However, at position two, 



c 



b^a 



> I 



will increase. Finally, a pole in position three near c , 



bs c 

-= <1 M is reduced while n is increased. Similar con- 

b^ a p 

elusions may be deduced when the poles are complex , 

The effect of a real zero at various positions is evident from 

Figure 4,4, As the zero approaches the point 1,0 the ratio becomes 
b c 

larger and increases while reducing n . A zero placed 

D a p 

very near the point 1,0 is not desirable since M increases tre- 
mendously, For most cases, it has been found that one zero placed 
between the origin and a position somewhat less than the real part of 
the predominate pole will result in a satisfactory transient response. 
This corresponds to using a phase lead network in conventional design 
methods. Other poles and zeros near the origin tend to cancel each 
other and have little effect on the time response. Thus, a pole- zero 
pattern, such as shown in Figure 4,5 consisting of a single real zero 
and a pair of predominate poles should serve as the starting point in 
the synthesis of transient response. 

The quantity n can be readily determined by considering the 
s 

second term in Equation 4, 10. Since the cosine term is always less 
than unity, the output will be certainly less than 5% of its final value 
value when the following equation is satisfied. 



n 

hj ' 



o 

• 

o 


2 K 






-PT^ 





(4.20) 



Taking the logarithm of both sides yields 
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Figure 4. 5. A basic pole-zero configuration for adjustment 
of transient response 
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1.3 - Log 



(4.21) 



n 

s 



Log 



K 



Q(^„) 

^o> 



z 



A response with no output at zero time requires the coefficient 
of the second term in Equation 4. 1 0 to equal the reciprocal of the 
cosine term with angle of 

[ ang Q(z^) - ang (z^ - 1) - ang P'(z^^) ] (4.22) 



For the pole-zero pattern of the type shown in Figure 4.5, the angles 
are those shown in Figure 4. 6 and generally the expression 4. 22 is 
about 180° and the coefficient term is nearly unity. In this case Equa- 
tion 4. 21 can be replaced by an approximated but very useful result 



n 

s 



1, 3 

log UqI 



(4.23) 



The solution of Equation 4.23 generally yields a noninteger value for 

n^ and the actual value is the upper integer value of the n^ obtained. 

The locus of z for constant n are circles on the z-plane with 
os ^ 

centers on the origin. The circles are shown in Figure 4. 7 and can 
serve as starting point for locating the predominate poles. 

The number of oscillations may be calculated approximately as 
follows 



or 



_ > settling time (4. 24) 

period of damped oscillation 

n 0 

N = - ? - ° (4.25) 

2 TT 
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z> plane 




Figure 4. 6. Various important angles 
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z- plane 




Figure 4. 7. 



z circles of constant n 
o s 
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For control systems with steady-state error to a unit step input, 
the constant term in Equation 4, 10 differs from unity by the extent 
of this error. The value of K may still be calculated since the amount 
of the error is given in the specifications. However, in most practical 
cases, the error is negligibly small and the same procedure previously 
discussed may be followed. 

The error constants of interest are the position, velocity and 
acceleration constants. These were developed in section 2,4 in terms 
of G(z) and repeated below for convenience, 

K* = lim G(z) 

P z-^l 

K* = 4^ lim (z - 1) G(z) 

V T ' 

z — >1 

K* = lim (z - 1)^ G(z) 

^ T'^ z -^1 

Introducing a change of variable 

6 = z - 1 (4. 26) 

The error coefficients become 

K* = lim G(0) (4.2 7 ) 

P e -^0 

1 lim 0 G(0) (4.28) 

e -^0 

1 2 

K* = — 5 — lim 0 G(0) (4. 29) 

^ T 0 -PO 

In terms of the variable 0 Equations 4. 2 and 4. 3 are 
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(4.30) 



w 

Q(6) = rr (6 + 1 - q, ) 

k=l 
m 

P(6) = XT (e + 1 - pj (4. 31) 

i=l ^ 

From Equation 4. 1 , the open loop transfer function of the sys- 

tem shown in Figure 4, 1 with H(s) equal to unity is 



w 



G(6) = 



K (e + 1 . 



m w 

7T (0 + 1 - p.) - K 7T (0 + 1 - qj^) 
i=l k=l 



(4. 32) 



Substitution of Equation 4, 32 into Equation 4, 27 and taking the limit 
yields 



w 








m 


w 


rr (1 . p ) . 

i=l ^ 


K TT (1 . q ) 
k=l ^ 



(4.33) 



Increasing the value of K will effect an increase of 
second approach would introduce a pole near the position 1, 0 and a 
zero between it and the origin. 

Incorporating the substitution 



a. 

1 




(4. 34) 



b 



k 




(4.35) 



into Equations 4.30 and 4,31 yields 
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w 



0(6) 


= TT 


(e t ^ ) 


(4.36) 




k=l 








m 


(e * ^ ) 




P(0) 


= TT 


(4.3 7) 




i=l 


1 





Introducing the change of variable and the substitution above in Equa- 
tion 4. 1 8 give s 

w 



K 



TT 

k=l 



7T 



i=i 



a. 

1 



(4.38) 



Thus it follows, for a system with unit feedback and no steady-state 
error to a position input, the closed- loop function becomes 



V^(0) 



w 

m 

TT (a. e + 1 ) 

i=l 



(4.39) 



The open- loop transfer function is then 

w 



TT (b, e + 1 ) 

G(6) = 

m w 

jj- (a. 0 + 1) - TT (b e + 1 ) 
i=l k=l 



(4.40) 



Expanding the above eqviation 



74 



I 



I 






/ tel t j 



G(0) 



V./ .2 V4/ ■ 

its / / /f- / 



W 4 / Vi/ 

+ ••' terrb ^ 

k^i 



d 



m vs/ 


e% 




^ 1 




^ \ 



^ ,2w f ^ £ 

ita/ 1 /i*/ 



-I- • • *f e"^ ^ ^i 
£-1 



nn 



/ 



/ 



Multiplying the above equation by and then setting Q to zero as 
indicated by Equation 4,28 yields the velocity error constant 



T K* 



(4.42) 



” m 

“ ■ J. 



W 



m 



If 



I ^ is zero the velocity error constant is infinite 



k=l i=l 

and the acceleration error constant becomes 



K* 



w m 

E 

fe^l i=l . 



(4.43) 



In terms of the z -plane poles and zeros 

1/T 



K* 



m 



y ^ 

^ i-p. 



w 

V 1 



i=i - Pi kti 



(4.44) 
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(4.45) 



The denominator of the above equation can be reduced (large 
by placing a zero near the point 1,0 . However, this will have the 
additional effect of increasing M and n^ , This will be illustrated 
by the following example. Consider a simple system with a configura- 
tion such as shown in Figure 4. 5 , 



K (z - q^) 



W(z) = 



(z - + j v^) (z - - j v^) 



(4. 46) 



Placing these poles and zeros of the closed-loop function into 
Equation 4,46 , placing the expression to zero, and solving for q^ 
which would yield an infinite velocity error constant gives 



q 



o 



2(1 . u^) 



(4.47) 



Assigning values of u^ = 0, 8, v^ = 0. 4 results in a q^ of 
0, 5 , This location for the zero will yield an infinite velocity constant 
but the overshoot, M , for this system is 0, 6 , If the zero had 
been placed at the origin the value of M would be a more reasonable 
0,39 but the value of T would be unity. Thus it appears that 

compensating with one zero does not give independent control over 
peak overshoot and velocity error constant. 
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An arrangement which does permit adjustment of 
appreciable effect on the transient response makes use of a pole and 
a zero placed very close together in the form of a dipole and near to 
the point 1,0 as shown in Figure 4, 8 . The zero at and the pole 
at p^ are placed sufficiently close together that the residue of the 
transfer function in the pole at p^ is made negligible. Although p^ 
and q. are nearly equal, there can be a significant contribution to 
the denominator in Equation 4,44 because the inverse values 

and - 5 — i are large numbers,. 

1 - q 1 - p 

Consideration of the Equation 4.45 shows that a dipole, but with 
the zero closer to the origin, may be used to raise the acceleration 
constant. Furthermore, additional zeros would be required for stabi- 
lization purposes giving rise to a typical pole-zero configuration shown 
in Figure 4.9 . 

4, 3 The Required Number of Poles and Ze ros 

The general location of the poles and zeros to achieve various 
requirements was shown in the previous section. However, at the start 
of the synthesis procedure it is necessary to know the minimum num- 
ber of poles and zeros that are needed by a given application. Although 
mathematically permissable, a practical compensator cannot possess 
poles at infinity. Consequently, the excess of poles over zeros for the 
compensated or desired open-loop function must be at least equal to 
the excess for the plant • Since the closed- loop function zeros 

are the same as those of the open- loop and the number of poles are the 
same for either open or closed loop function, it can be said that 
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Figure 4.8. Pole-zero configuration for large velocity error 
constant 




Figure 4.9. Pole- zero configuration for large acceleration 

error constant 
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(number of poles - number of zeros) 



W(z) ^ 

(number of poles - number of zeros) ^ 

^p(z) 



(4.48) 



It would be of considerable advantage to avoid the lengthy pro- 
cess of taking the z-transform from the L*aplace transform in order 
to determine the number of zeros of the plant. In general this can be 
done by inspection as will be illustrated by considering the open- loop 
systems shown in Figure 4, 10, When the continuous system as 
shown in Figure 4, 10a has two poles in excess of the zeros, output 
will be available at the second sampling, i, e, , n = 1, Since the value 
of the output can be determined at the sampling instant by dividing the 
denominator of G(z) into the numerator, an output at this sampling 
instant requires that the number of z-plane zero be only one less than 
the number of z-plane poles. The number of zeros equals the number 
of poles when the Laplace transform poles exceed the zeros by only 
one. Also, in these cases, a zero is found at the origin. If a hold 
network (1) is employed, as shown in Figure 4, 10b, the same pro- 
cess of reasoning indicates that the number of poles exceeds that of 
zeros by one; however, a zero need not necessarily be at the origin, 

A plant including a delay, such as Figure 4, 10c , of less than one 
sampling period will give rise to the same z*wplane pole-zero differ- 
ence as discussed for Figure 4, 10a except the zero will not gener- 
ally be found at the origin. For every complete sampling period that 
the delay involves, poles are placed at the origin without any additional 
zeros. The overall transfer function for Figure 4, lOd is the product 
Gpi (z) (z) and the difference between the poles and the zeros 
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(a) The plant is a continuous network 



t(S) 








CCS) , 



(b) The plant includes a hold network 



^(S) ^ g(z-) 




C(g) 

C(^) 



(c) The plant contains a pure delay 



CCi) 




CCS) 



(d) Parts of the plant are separated by samplers 



Figure 4. 10. Various types of plents 
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of each function is additive in the final transfer function. If the poles 
and zeros differ by one, then the over-all transfer function will have 
an excess of two poles over the zeros. This is in accordance with phy- 
sical expectation since the second network does not receive an input 
until the second sampling instant. 



4.4 Open- Loop Root Locus 

The method of this dissertation assumes the closed-loop function 
to be known leaving the open- loop function to be determined. For a 
system such as Figure 4. 1 the open-loop function G(z) and the 
closed- loop function W(z) are related by 

= -TT^ 

Solving for G(z) gives 



G(z) 



W(z) 

1 -'W(z) 



(4. 50) 



Hence, the poles or roots of G(z) 

1 - W(z) = 



are the solution of the equation 
0 (4.51) 



Both functions are rational polynomials in z so that 



W(z) 



K Q(z) _ 

■PliT ■ 



K TT (z - q ) 
k=l 



7T (z - p.) 

i=l ^ 



r 



G(z) 



k A(z) 
B(z) 




(z - Cj^) 



m 



TT 

i=l 



(z 



Pi) 



(4. 52) 



(4. 53) 
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Equation 4, 50 may be also written as 



kA(z) _ K Q(z) 

-bTTT = P (z) - k Q (z) 



(4. 54) 



The gain constant are the same, the zeros of the open- loop function 
are identical with the zeros of the closed-loop function and the poles 
of the open-loop function are the value of z for which 



W(z) = e-* ^ ^ ^ (4. 55) 

The locus of roots consist of all points on the z- plane for which the 
plane of the closed-loop is 2 h tt (h may be any integer or zero) 
that is, 

ang [ W(z) ] = 0 + 2 h tr (4. 56) 



In terms of Equation 4, 52 the phase of vV(z) is given by 
w m 

0 = ^ ang (z - qj^) - T ang (z - p.) (4. 57) 

k=l i=l 



The precise position of the poles on the root locus requires that 



I W(z) I 



(4.58) 



or that 



w 



K TT 

k=l 



TT 



i=l 





(4. 59) 



The first step in applying this graphical root finding method is 
to determine the locus of roots. An approximate root locus can 
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generally be found very quickly and simply by applying some simple 
rules concerning the geometry of the zero phase lines in the z-plane. 

7 

After the approximate locus has been obtained the exact location of the 
roots can be determined by use of Equation 4. 57 and 4, 59* • A typi- 
cal open-loop root locus is shown in Figure 4, 11 

The rules and their proofs for the approximate construction of 
the open-loop root locus are now presented. 

Rule I , ”The loci are continuous curves starting at a pole for 
zero gain and terminating at a zero at infinite gain, 

Rule 2, **The loci exist on a part of the real axis where an even 
number of poles plus zeros or none are found to the right, ” 

This may be easily verified by consideration of Equation 4, 55 
and Figure 4,11, 

Rule 3 , If the number of poles, w, exceeds the number of zeros, 
m , there are w - m branches which go to infinity. Further, the 
directions of the infinite asymptotes with respect to the real axes are 



2 h tr 
w - m 



(4,60) 



where n = 0 , 1,2, 3 

This rule is proven by reference to Figure 4, 12a and b. The 
scale of the axes in Figure 4, 1 2b is increased over that of Figure 4, 1 2a 
reducing the location of the poles and zeros of W(z) to almost a point. 
As point P moves further away from the origin the angles subtended 
by the zeros or poles are essentially the same. The sign of the angles 
from the zeros is apposite to the sign from the poles. Since the angles 



*For this operation a device called the Spirule has been developed by 

Spirule Company, 9728 El Venado, Whittier, Calif, 
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Figure 4. 11. A typical open- loop root locus 
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(a) 




(b) 

Figure 4. 12. Proof of rule three 
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from the excess poles are all the same, the angle of each pole is 



if the sum is equal to 2 ir . The following table will 

w - m ^ ® 



prove useful. 


TABLE 4. 1 

Direction of Asymptotes 


w - m 


Directions 


1 


0 


2 


1 

00 

o 

o 

o 


3 


-120° 0 +120° 


0 

o 

1 


0 +90° + 180° 



Rule 4. 


’’All asymptotes intersect at a point u^ on the real 


axis. " 


m w 


where 


^ p. - Z-i q, 

i=i ^ k=i 

u = (4.61) 

w - m 



The proof of this rule can be demonstrated by consideration of 
the transfer function of W(z) written in the form 



W(z) = 


m , m-1 , 

K z + a^z + . . . + a^ (4. 62) 

W . u W-1 

z + bj^z + . . . + 



Dividing the numerator into the denominator to give 



W(z) = 


K 

w-m . w-m-1 , (4.63) 

z + (bj - a j) z + . . . 
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The locus is described by equating the denonminator to - K, that is 



w-m 

z 



(bj- 






w-m- 1 



+ . 



- K 



(4. 64) 



By the theory of equations the sum of the roots of Equation 4, 64 is 
bj^ - a^ . When the magnitude of z is large Equation 4. 64 behaves 
as a polynomial of w-m degree with the sum of the roots given by 
b^ - a^ being a constant which is independent of gain as long as 
w-m is greater than two. This is no limitation as an inspection of 
Table 4, 1 will show. With independent of gain, the sum of 

the root remains at one point for any value of gain. Hence, the asymp- 
totic lines radiate from the point on the real axis given by 




w-m 



(4.65) 



Rule 5. "The point u^ at which the locus breaks away from 
the real axis is found by proper equating of the sum of reciprocal of 
the distances between poles and zeros and the point u^ as explained 
by the following: " 

The appreciation of this rule and its proof is demonstrated by 
Figure 4. 13. A trial point is assumed at u and moved a small 
vertical distance A v off the real axis. To be a point on the root 
locus this point must satisfy the phase angle criterian. That is 



- e, + 0- + e, +0. = 2 h ir 

12 3 4 



(4.66) 



All the angles are small and can be replaced by their tangents, hence 
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Figure 4. 13. Proof of Rule Five 
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A V 



+ 180° + 180° = 360° (4.67) 

u. - u ' ' 



which becomes 




u 



3 



I 



e 




I 



(4.68) 



The solution of Equation 4.68 for the value of u^ is generally 



accomplished by trial and error. For the case where complex poles 
or zeros exist, the expression needed to determine the break away 
point is more complicated. However, it has been shown (21) that unless 
the complex poles or zeros are close to the point of emergence their 
contribution is slight. 

Rule 6. ’’The angle at which a locus emanates from a complex 
pole or terminates upon a complex zero is found by summing the angles 
to all the other poles and zeros. When emanating from a pole the 
angles from the poles are taken negative. When terminating on a zero 
the angles from zeros are taken negative.” 

Consider a test point P shown in Figure 4. 14 as on the locus 



only a very small distance removed from the complex pole. Then the 
angle can be considered the angle of departure of the locus from 



the pole. Since P^ is on the locus the sum of angles must apply. 



64 + - (63 + e +65 + 0^) = 0 



(4. 69) 



solving for 



«6 = ( - («3 + «4 + > 



( 4 . 70 ) 



As point approaches the pole, the angles to the right of the 
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I 

I 



above equation become simply the angle subtended between the various 
poles and zeros and the complex pole in question. Obviously, when 
the angle of arrival at a zero is determined a change of sign is required. 

4. 5 Closed- Loop Root Locus 

Application of the closed- loop root-locus assumes that the poles 
and zeros of the open- loop function are known and that the poles of the 
closed- loop function are to be determined. For a system such as Fig- 
ure 4, 1 the open-loop and closed- loop function are related by 

= 1 t 

Hence the poles or roots of W(z) are the solution of the equation 

1 +G(z) = 0 (4.72) 

As discussed in the previous Section, both functions are rational poly- 
nomials in z . Further, the zeros and gain factors of both functions 
are identical. Thus, the poles of the closed-loop function are the 
values of z for which 

r-t \ j h ir 

G(z) = e^' 

h = 1, 3, 5 , (4.73) 

The precise position of the poles or the root locus requires that 

I G(z) 1 = 1 (4. 74) 

As in the case of the open-loop root locus, the first step in apply- 
ing this graphical root finding method is to determine the locus of roots. 
The approximate root locus are easily found by application of the rules 
which have been presented (6) (15) in the literature. A typical closed- 
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loop root locus is shown in Figure 4, 15, 



4, 6 Additional Design Aids 

It has been tacitly assumed, before the techniques of this disser- 
tation are employed, the optimum system performance by simple gain 
adjustments was found short of the required specifications. 

This preliminary gain adjustment can be aided by Figures 4, 16 
and 4,17, These Figures relate the transient response quantities of 
M and n^ for second order systems with no position error to the 
location of the closed- loop complex poles. Figure 4, 16 applies to a 
system with sampler only, while Figure 4, 17 applies to a system incor- 
porating a zero order hold network (1) (12) , Only the first quadrant 
is shown since complex pole- zero patterns are symmetrical about the 
real axis and poles in the second quadrant seldom give rise to a satis- 
factory response. The figures are entered on the circle correspond- 
ing to the product of s -plane open- loop pole and sampling period. These 
circles are actually the closed-loop root locus of the complex poles 
and proceeding counter clockwise corresponds to increased gain. On 

the circles are indicated the values of ,M and n that a second 

P 

order system with a pole at that location will produce. For complex 
systems, the figures can be used on the basis of dominance to deter- 
mine if the various values of M and n^ will satisfy the specifica- 
tions with sufficient margin to provide for the deleterious effect of the 
neglected poles and zeros. If this appears possible, a more exact 
study can be made using the closed-loop root locus. The figures may 
also be helpful for making the first selection of the predominate poles 
in synthesis process. 
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z -plane 




Figure 4, 15. A typical closed- loop root locus 
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Figure 4. 16, Values of M and for second order systems 




Figure 4, 17, Values of M and n for second order system 
with hold P 
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The dipole required to raise the steady- state error constant may 
require such close spacing as to render praphical methods impractical. 
Applying the procedures of this thesis analytically would require the 
finding of roots of polynomials. Since the approximate values of the 
roots are already known, this is not as difficult a task as may first 
appear. However, an alternate method avoiding these difficulties works 
directly with the open- loop function. A pole and a zero are introduced 
into the open- loop function of proper relative values to give the desired 
increase of steady-state error constant as indicated by Equations 2.64, 
2, 65 and 2. 66, The expression for the output transform can be found 
as 



C(z) 



z k A(z) 

z~- r ^ k A(z) + 6(z) 



(4, 75) 



where A(z) and B(z) are the numerator and denominator of the 
open- loop function while k is the gain. The actual output at the sam- 
pling instants are readily found by division. 

Unfortunately, this method leaves the designer in doubt as to 
whether the final output time response will be satisfactory. Actually 
the above procedure is not so involved that trial and error is out of 
order, Furthe r, as an aid in the initial selection it can be shown that 
the closer the loop pole of the dipole is to unity the smaller will be the 
residue of the output transfer function in the pole. 

This effect of open- loop pole location on transient response when 
employed to increase the steady-state velocity error coefficient will 
be illustrated as follows. For simplification let the original value of 

and the sampling period be both equal to unity. Let it be required 
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that the value of K*^ be raised by a factor of five. If the open-loop 
pole and zero are designated by a and q and the closed- loop pole by 
p then it follows directly from Equations 2.65 and 4.44 that 



5 = 



1 - q 

1 - a 



(4. 76) 



and 



5 



1 + 



~r 

i-p 



~T" 

1-q 



(4. 77) 



Solving Equation 4. 77 for p yields 



P 



9 q - 4 
4 q + 1 



(4. 78) 



The introduction of a dipole into the closed- loop function affects 
the response in two ways. 

(1) The original terms are modified, but slightly since, as a 
result of the fraction (z-q)/(z-p) which is essentially unity at the 
original poles, the residues at these poles are almost unchanged. 



(2) The term (p 1 ^ P ' ^(p) added . 

Consideration of this term will show that with movements of the 
dipole near the point 1.0 the coefficient of this term will change by 
a factor almost directly equal to (q-p)/(p-l) . Figure 4. 18, 

employing the previous expressions, gives the value of this ratio as 
the open- loop pole is moved closer to unity in value. 
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m 



Relative magnitude of residue 







Open- loop pole a 

Figure 4. 18. The effect of the position of the open- loop pole 
on the output 
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CHAPTER V 



APPLICATION OF SYNTHESIS PROCEDURE 
5, 1 Introduction 

Chapters III and IV are concerned chiefly with presenting the 
synthesis philosophy and the techniques necessary to implement the 
philosophy. In this chapter the actual process of synthesis is illus- 
trated by a variety of different specifications. Also discussed in 
detail is how the methods of this dissertation may be applied to sys- 
tems which do not have a sampler in the error channel. The chapter 
is concluded by a discussion on synthesis of systems with multiple 
inputs. 



5.2 Complete Example of Synthesis Procedure 

This section presents an example of synthesis using the method 
of this dissertation, an example purposely chosen to illustrate the 
basic elements of the procedure. Synthesis is initiated by determin- 
ation of the specifications which include the fixed part of the system 
and the time domain performance characteristics. The specifications 
are as follows: 

(1) The plant has a sampler in the error channel operating with 
period of 0, 1 second and a transfer function given by: 



Gp(s) 



K (s 19) 
s(s+3. 7)(s+8)(s+25) 



(2) The time domain specifications: 



(5. 1) 
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Approximating the plant as a second order system and employ- 
ing Figure 4. 1 5 as explained in Section 4. 6 shows that specifications 
cannot be satisfied by adjustment of gain and a compensator will be 
necessary. Consequently, the synthesis procedure outlined in Figure 
3. 3 will be required. 

Two items which must be decided before selection of the closed- 
loop pole-zero configuration are the difference between the number of 
z-plane poles and zeros of the plant and how many of the poles and zeros 
of the plant are to be retained in the desired transfer function. For 
the first item it is seen that the plant has one sampler, no delays and 
a numerator three degrees less than the denominator, hence the z- 
transfer function of the plant will have one more pole than zero. For 
the second item, the retention of plant poles and zeros is largely tied 
up with consideration of the physical equipment. In this case an infin- 
ite value for requires the pole at the origin to be retained and 

arbitrarily the pole at 3. 7 will also be retained. 

With these constraints the simplest pole-zero configuration can 
be a pair of complex poles and one real zero. Usirt^- Figure 4. 1 5 
again as a guide the first choice of pole-zero selection is made with 
poles at 0. 5 j 0.45 and a zero at 0. 3. Using the correlation 
theorems 

M = . 21 



n = 8 

s 

K* = 13 

V 



The transient response requirements are met with sufficient mar- 
gin to permit realization of a higher velocity error constant by intro- 
duction of a dipole. First, however, the open- loop poles must be ascer- 
tained. This is readily done by the open-loop root locus shown in Fig- 
ure 5. 1. The gain is fixed by including one pole at point 1. 0 and then 
the remaining pole is found at 0. 64. Also shown in this figure are the 
construction lines used for measurements and angles to determine n 

P 

and M . n is estimated by reference to Figure 4. 7 . 
s 

The s -plane pole of 3. 7 which is to be included in the final open- 
loop function corresponds to a z-plane pole of 0. 691. Since the first 
try has produced a pole not too far removed from this location, a 
closed- loop root locus is now formed based on the zero at 0. 3, and 
poles at 1 and 0.691 . This locus is shown in Figure 5.2 . Choos- 
ing on the locus poles at 0. 55 j 0. 46 gives the following time 
% 

specifications 

M = .28 

n = 3 

P 

= 9 

= 13 



With the gain determined, the open-loop function is 



G(z) 



. 591 (z - .3) 
(z - I ){z - . 691 ) 



(5.2) 



In accordance with the discussion in Section 4. 6 an open-loop 
pole is added at 0.99 and a zero at 0.95. The value of K* 

V 



100 



z-plane 

X closed-loop pole 




Figure 5. 1, First selection for pole-zero configuration and 
open-loop root locus. 
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Figure 5, 2, Second selection for pole-zero configuration and 
closed-loop locus. 
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then given by 



K. . -591 . 65. (5.3) 

(1 - .99)(1 - .691) 

Further, realization of the open-loop function as a continuous function 
without time delays requires a zero at the origin. To retain the dif- 
ference of one between the number of poles and zeros an additional 
pole is added at 0.001. This pole was purposely placed very near 
the origin to cancel the effect that the zero may have on the transient 
response. The desired open-loop system is thus described by the 
transfer function 



G(z) = K ( 5 . 4 ) 

B(z) (z-l)(z-. 99)(z-.691)(z-.001) 

The inclusion of the additional poles and zeros will result in a 
shift of the complex poles and a change in the transient response. The 
new pole positions could be determined and the correlation theorems 
employed once again. But a simpler method would be to form the out- 
put z-transform function C(z) to a unit step directly from G(z) as 

C(z) = — X KjMzj ^3 

z - 1 K A(z) - B(z) 

By simple (division the output at the sampling instant is in this case 

C(t) = . 591 u^(t-. 1) + 1. 087 u^(t-.2) + 1. 33 u^(t-. 3) (5.6) 

+ 1. 34 u^(t-. 4) + 1.212 u^(t-. 5) + 1.065 u^(t-. 6) 

+ 0. 966 u^(t-. 7) + 0. 935 u^(t-. 8) + . 955 u^(t-.9) 

+ . 995 u^(t-l.O) 

1 03 



} 









The transient specifications are then seen to be 

M = 1.34 

n =4 

P 

n =9 

s 

Converting the z-transfer function to its equivalent Laplace trans- 
form with poles in the primary strip is the next step. This involves 
taking a partial-fraction expansion of 

G(z) 

z 

i 

and multiplying through by z , converting term by term to Laplace 
transform with the aid of Table 2. 2 , and, lastly, collecting the terms 

over a common denominator. These three steps take the form shown 
below. 



G(z) 



6. 698z 
z - 1 



5. 514z 
z - .99 



. 9382z . 246z 

z - . 691 z - . 001 



G(s) 



6. 698 
s 



5. 514 _ . 938 _ . 246 

s + .l s + 3. 7 s-6. 9 



20. 95(s -I- . 512)(s -I- 15.94) 
s(s + . l)(s + 3. 7)(s + 69) 



The compensation is now determined by dividing the derived 
transfer function. Equation 5.9, by the plant Equation 5. 1, yielding 



2 0. 95(s + . 512)(s+8)(s-H5. 94)(s+25) 
(s + . l)(s + 19)(s + 69) 



(5. 10) 
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The network is not realizable in its present form but the addition 
of a fourth pole far out on the negative axis, has a negligible effect on 
the transient response, and yields a realizable compensation function 



ZO. 95 (s+, 5IZ)(s+8)(s-H5, 94)(s+Z5) 
(s + . l)(s + 19)(s + 69)(s + 100) 



This function is realizable in terms of resistances and capacitors by 
standard network synthesis procedures (IZ) (Z3) . 

The preceding paragraphs have attempted to illustrate the basic 
synthesis procedures. However, experience, personal preferences 
I and considerations of particular systems may indicate deviations from 
the initial steps outlined above. As for example, if the entire plant 
must be retained as part of the final desired open-loop function it may 
be expedient to begin the design on that basis. This would correspond 
to entering the Flow chart (Figure 3. 3) initially from the plant posi- 
tion, This versatility of design procedure is the strong point of this 
synthesis method. This versatility is furnished by the two root locus 
techniques which relate the movements of poles in one function on the 
other and the correlation theorems which permit interpretations of the 
pole-zero configuration in terms of time response. 



5,3 Further Examples of Syste m Design 

In this section three systems are considered, each with the same 
continuous fixed components involving one integration and an s-plane 
pole at 0,431, The first system includes only the required sampler, 
the second employs a sampler and holding network, while the third is 
characterized by a pure time delay of one half second, giving rise to 
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the following three transfer functions. 






K 



s(s + . 431) 



(5. 12) 



^2p(^) = 



1 - e 



-s 



K 



s(s + ,431) 



(5. 13) 



G3p(s) = 



K 



s(s + . 431) 



X e 



- .5 s 



(5. 14) 



In each case the sampling period is one second. 

Besides requiring that the above fixed elements be included in 
the final open- loop transfer function each system must also satisfy 
the time specifications. 



M 
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n 
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K 

V 




1 



With all systems required to meet the same time specifications 
at the sampling instants and having the same principle s -plane pole, 
the same closed-loop z-transform should be satisfactory for all three 
systems. A pole-zero configuration with one pair of poles at 0. 5 it 
jO. 55 and a zero at 0. 15 meets the above specifications w^th 



M 


= 


.33 


n 




2. 56 


n 


— 


3 


P 






K 


= 


1.3 



V 



Forming the output transform to a step function and dividing out 
gives the output time sequence as 
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the following three transfer functions. 



G 
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(s) 



K 

s{s + . 431) 



(5. 12) 
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s{s + .431) 



{5. 13) 
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3p 



(s) 



K 
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s(s + . 431 ) 



-.5 s 
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(5. 14) 



In each case the sampling period is one second. 

Besides requiring that the above fixed elements be included in 
the final open-loop transfer function each system must also satisfy 
the time specifications. 
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With all systems required to meet the same time specifications 
at the sampling instants and having the same principle s-plane pole, 
the same closed-loop z-transform should be satisfactory for all three 
systems. A pole-zero configuration with one pair of poles at 0. 5 i 
jO. 55 and a zero at 0. 1 5 meets the above specifications w^th 
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Forming the output transform to a step function and dividing out 
gives the output time sequence as 



106 



c(t) 



= 0. 65 u (t- 1) + 1. 202 u (t-2) + 1. 408 u (t-3) 

o o o 



+ 1.306u (t-4)+. . . 
o 



(5. 15) 



The reason for the difference between the predicted value of M and 
its actual value is due to the noninteger value for n as explained in 
Section 4. Z . 

The closed-loop transfer function has a corresponding open- 
loop transfer function 



GU) , --.Asu --^1> 

(z - 1 )(z - . 65) 



(5. 16) 



It should be noted that the plant pole has been accommodated 

since 

" -^31 

e = . 65 

Although to this point the procedure has been the same for all 
three cases, the conversion to the l^place transfer function, the next 
step in deriving the compensating network, is different for each case 
and will be illustrated in turn. 

(1) For the first system the plant is a continuous network whose 
z-transform is characterized by a zero at the origin. Adding this zero 
to Equation 5. 15 and a pole 0. 001 to maintain the proper difference 
between poles and zeros, the resulting transfer function becomes 



Gid (z) = 



. 65 z(z - . 15) 



(z - l)(z - ,65){z - .001) 



(5. 17) 



Employing the methods described in the previous section yields 
the equivalent Eaplace transfer function 
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I 



1. 6(s + 2. 82) 
s(s + . 431)(s + 6. 9) 



(5. 18) 



G^^(s) = 



the transfer function of the compensator is then 



Gd(s) 



0.408(s + 2. 82) 
(s + 6.9) 



(5. 19) 



(2) For a system including a zero order hold network the trans- 
fer function of the open-loop system is the general form. 



— K G G (z) (5.20) 

s o 
z 

The (z - 1)/ z term corresponds to the transfer function of 
- s T 

1 - e , while the transfer function G^(s) of the pure integra- 

tion supplied by the hold network and the transfer function G^(s) of 
the continuous portion of the system result in the z-transform 
G^ ^*20 form of Equation 5. 20 yields 

G (s) = -- - X - (5.21) 

z (z - 1) (z - . 65) 

and then 

G G (z) = ~ ■ (5.22) 

(-1)^z-.65) 

The conversion of G^ ^ZO ^ corresponding Laplace trans- 

form is similar to the previous example yielding 
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(5.23) 



G 

s 



(s) G 



20 



(s) 



. 449(s -I- 1.514) 
(s + .431) 



hence 



8. 98(s + 1. 514) 
s(s + . 43 1 )(s + 20) 



(5.24) 



Where the pole at 20. 0 has been added to enable physical reali- 
zation of the compensator 



G 



2c 



(s) 



0. 0757(s + 1. 514) 
(s + 20) 



(5.25) 



(3) The third system includes a pure time delay of one-half the 
sampling period. Because of this delay a zero at the origin is not 
required in the transfer function. Proceeding directly to the partial 
fraction expansion of Equation 5. 16 gives 

G (z) = - + — (5.26) 

z - 1 z - , 65 

The conversion of each term on the right hand side of the above 
equation to the s-domain must be done through use of Table 2,3 , 

Modified z- transforms , with m equal to one-half, giving rise to 



^3d _ 1. 58 ^ 1. 15 

-.5s s S+.431 

e 



(5.27) 



Combining terms under a common denominator, adding a pole at twenty 
and making the necessary gain adjustment gives 
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8. 58(s + 1.52) 
s{s + . 43 l)(s + 20) 



-.5s 



(5.28) 



G 



3d 



(s) 



e 



the required compensator is 



G 



3c 



(s) 



0. 075(s + 1. 52) 
(s + 20) 



(5.29) 



A second series of three applications evolves when the value for 
is increased to six while holding all other specifications the same. 
The open- loop z-transfer function satisfying these conditions is 



G ^ (z) 



. 65(z - ,92)(z - . 15) 
(z - l)(z - . 98)(z - .65) 



(5.30) 



Employing previously discussed methods of this chapter the 
required open-loop systems are described by the transfer functions. 



1. 574(s + . 083)(s ^ 2 . 88) (5.31) 

s{s + .02){s + .431)(s + 6,9) 



^2d<^) 



1 - e~^ ^ 9. 477(s + . 8)(s + 1. 43) 

s s(s + . 02)(s + . 431)(s + 20) 



(5. 32) 



G 



2d 



(s) 



8. 45(s + . 083)(s + 1. 572) 
s(s + . 02)(s + .431)(s + 20) 



(5.33) 



The output time sequence to a unit step input for the above sys- 
tem is found to be 



no 









I 



\ 

I.I 



(5.34) 



c(t) = 0. 65u^(t-l) + l,24u^(t-2) + 1.^8u^(t-3) 

+ 1. 37 u (t-4) + . . . 

0 

Each of the six systems synthesized in this section were instru- 
mented on an analog computer. The computer outputs are given in 
Figures 5.3 through 5,8, It is to be noted that the output at sampling 
instant checks very closely with the calculated output. iSble 5. 1 com- 
pares the calculated and measured values for 

Liberty was taken in Figures 5, 5 and 5, 8 where the curves are 
actually the input to the delay network displaced by half a sampling 
period. The imperfection of the computer time delay required this 
modification of the actual data. Justification for this procedure as 
well as detailed description of the computer study is given in Chapter 
VI. 

5,4 Sampled Data System without a Sampler in the Error Channel 
The sampled data systems to be discussed in this section are 
characterized by a lack of sampling in the error channel. Figure 5.9 
is the block diagram of such a system. Reference to Table 2,2 gives 
as an expression for the output 

R G, (z) (z) 

C(z) = 1 ^ (5,35) 

1 + H Gj G^ (z) 

It is of particular interest that the transform of the input signal does 
not appear as a multiplicative factor in the numerator of the expres- 
sion for C(z) . Consequently, it is not possible, as was done pre- 
viously, to divide out the transform of the input signal and form a 

closed-loop transfer function independent of the input, 
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Figure 5. 4. Response of system with hold to step input 
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Figure 5, 6. Response of system with sampler only to step input 
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Response of system with delay to step input. 



Table 5. 1, Comparison of calculated and measured 
velocity error constants 



System 


Calculated 


Measured 


sample r 


1. 58 


1. 68 


hold 


1. 58 


1. 64 


delay 


1. 58 


1. 55 


sample r 


6.2 


5.4 


hold 


6.2 


5. 7 


delay 


6.2 


6.0 
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Figure 5. 9. 



A sample-data control system 



By the expediency of dividing both sides of Equation 5. 3 5 by a 
z-transform of the input a closed-loop transfer function can be defined 
as 

P R G, (z) (z) / R(z) 

-^(z) = i (5.36) 

I + H G^ G^ (z) 

It should be noted that the transfer function applies only for a specific 
input and for this reason will be referred to as a limited closed-loop 
transfer function with designation (z) . 

In general, limited closed-loop transfer functions do not have the 
same zeros as the open-loop transfer functions even though the system 
provides a unity feedback path. Because of this characteristic it is 
not possible to employ the open- loop root locus for finding the poles of 
the open-loop functions nor are the expressions for the steady state 
error constants in terms of closed-loop poles and zeros applicable. 
However, the response to any particular input can be evaluated dir- 
ectly in the manner employed in Section 4. Z for a unit step and the 
response shaped as desired by proper selection of the closed-loop poles 
and zeros. The steady state error can be found by direct application 
of the final value theorem. Equation Z.55. For systems with zero 
steady state error to a step input. Equations 4. 14 and 4. I 9 are 
directly applicable for synthesis to a prescribed transient response. 

With the limited closed-loop transfer function determined, the 
desired open- loop transfer function may be found by a method such as 
the following. With a step input and setting H(s) to unity the remain- 
ing expressions defined in Equation 5.36 are polynomials in z with 
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expressions 



R(z) 



z - 1 



(5. 37) 



R (z) = 



k , a , (z) 
rl rl 

(z-1) bj (z) 



(5.38) 



(z) 



^2 ^2 

(z) 



(5.39) 



Gj G^ (z) = 



bj (z) b^ (z) 



(5.40) 



(z) 



K Q (z) 

P (z) 



(5.41) 



In Equations 5,38 and 5,40 , which are z-transforms of the product 
of two Laplace transforms, the convenient property that s-plane poles 
have a one to one image on the z-plane is utilized by factoring the 
denominator into the parts corresponding to each Laplace transform. 
Substitution of the above expression into Equation 5,36 yields 



K 



Q(z) 

T^) 



k , a , (z) k~a.,{z) 
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(5,42) 



Generally, either Gj^ (s) or G^ (s) are known in part or 
whole. Assuming that G^ (s) is fixed, then the zeros of G^ (s) 
are fixed by Equation 5,42 , after the closed- loop zrros and gain con- 
stant are selected. Consequently, the factors of b^ (z) must be 



chosen to satisfy Equation 5. 39 and also combine with the remaining 
open-loop poles and the open-loop zeros to form the poles of the closed- 
loop function with the exception of the pole at the origin introduced by 
the transform of the input, 

A first selection of the poles is aided by the consideration of the 
physical system and the type of compensation to be employed. A phase 
lead network requires a pole near the origin, A dipole for reduction 
of the steady state error requires a pole near the point 1,0 , The 
total number of open-loop poles is one less than the number of closed- 
loop poles. The expression for A^(z) can be solved for by equating 
the numerators of each side of Equation 5, 42 , A first trial expres- 
sion for (z) is then known and function (s) can be deter- 

mined by methods illustrated in the previous sections of this chapter. 
Then the expres sion for Gj^ G^ (z) may be formed and by root locus 
technique the closed-loop poles found. In general, these poles will 
not be those previously selected. However, the resulting poles may 
very well satisfy the specifications thereby completing the design. 

If this is not the case the process must be iterated. The precise 
method for determining G^ (s) may vary with different system 
arrangements, but with experience systems of this sort can be readily 
handled by application of the above principles, 

5, 5 Systems with Multiple Inputs 

Because of the particular application or because of disturbances, 
all feedback control systems have multiple inputs. For example. Fig- 
ure 5, 10 shows the block diagram of a sampled data servo system in 
which the output is corrupted by a spurious signal u(t) entering the 
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Figure 5, 10, Tandem single- loop system with secondary 
input 
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system at point A , Thus a feedback control system is characterized 
by n closed-loop transfer functions, the customary or primary trans- 
fer function relating output with the reference input and a number of 
secondary transfer functions for each secondary input to reference out- 
put. Whatever the source of the secondary input or wherever the 
point of entry to the system, satisfactory design involves the realiza- 
tion of suitable dynamic characteristics for each of the possible trans- 
missions. 

In Figure 5, 10 (s) and (s) may represent transfer 

functions of fixed parts of the system while G^ (s) is the transfer 
function of the compensator. For a signal u(t) with transform U(z) 
the expression for the transform of the output is 

U G (z) 

C(z) = ^ (5.43) 

1 + ^2 ^3 

All factors on the right side of the above Equation have been 
fixed by synthesis of the primary transfer function given in the speci- 
fications, No provision is available for modifying C(z) if the results 
of Equation 5,43 are unsatisfactory. More latitude is offered if the 
compensator transfer function G^ (z) is split into tandem and shunt 
sections. This arrangement is illustrated by the system in Figure 
5,11a , Figure 5,11b shows the same system reduced to a single 
loop. The transfer functions of the blocks are related by the equations 

G^(s) = Gj (s) (s) (5.44) 
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(a) System with two input employing feedback and tandem compensator 




(b) Tandem single- loop system 



Figure 5,11. Two equivalent systems 
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(5.45) 



G^(s) 



G 3 (s) 

1 + (s) G3 (s) G^ (s) 



The closed-loop z-transform transfer function for the system in Fig- 
ure 5.11b is given by 




M 

1 + (z) 



(5.46) 



Figure 5. 12a shows the system of the previous figure redrawn 
and the reference input omitted to emphasize the transfer character- 
istics of the system between secondary input and reference output. By 
block diagram manipulation the system can be reduced to that shown 
in Figure 5.12b where the transfer functions G^ (s) and G^ (s) 
are the same as defined by Equations 5.44 and 5.45 . The transfer 
function of the output is given by the expression 



C(z) 



U G 



b 



(z) 



1 + G G, (z) 
a b ' ' 



(5.47) 



Using the methods of the previous section the limited closed- 
loop transfer function for the secondary input is obtained as follows 



C 

U 



(z) 



U G^ (z) / U(z) 



1 + G G, (z) 
a b ' ' 



(5.48) 



Noting that the denominators of Equations 5.46 and 5,48 have a com- 
mon term permits as a more convenient expression for 
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(a) 




Figure 5. 12, Two equivalent systems 
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u 



(z) 



U (z) 



U(z) G^ (z) 



R 



(z) 



(5.49) 



Equation 5.49 places in evidence the pos sibility of satisfying 
independently specifications for both the primary and secondary trans- 



mission characteristics. However, U(z) and R(z) must be the 
same function if the equation is to be meaningful. Poles of 



U 



(z) 



and 



R 



(z) are identical with the exception of an additional pole 



at the origin of the expression for 



— (z) . Consequently, the 

C 



design freedom for meeting specification on (z) is limited to 

choice of gain constant and zeros of U G^ (z) . 

The choice of zeros and gain constant for meeting specifications 
are aided by the relationships given in Section 4. 2 . With the design 



of the system for 
C 



(z) already completed and with the zeros 



of 



U 



(z) selected Equation 5. 49 may be solved for the remain- 



ing unknown U G^ (z) , Hence 



u (z) = 



_C_ 

u 



(z) 



R 

C 



(z) 



U(z) G^ (z) 



(5. 50) 



By procedures described in detail earlier in this chapter the transfer 
function for G^ (s) maybe evaluated. The transfer function of the 
feedback compensator G^ (s) is then given by 



(s) = 



G 3 (s) 

Gb («) 



G 3 (s) G^ (s) 



(5. 51) 
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It is apparent that the poles of (s) and (s) , which are 

also poles of the open- loop function, are zeros of G^ (s) . Conse- 

quently G^ (s) will in general require one or more poles at infinity. 
Hence, in practice G^ (s) can only be realized approximately by R- 
C circuits and high gain amplifiers. 

The expression for the transfer function Gj^ (s) , given in 

terms of the previously determined desired open-loop transfer func- 
tion, Gj (s) , is 

a 

G , (s) 

G, (s) = (5.52) 

^ G^ (s) G^^ (s) 



5. 6 Suggestions for Future Work 

The field of sampled-data theory and application is still in its 
infancy. Many problems remain to be studied. One field that should 
prove very useful is the use of sampled-data theory in the analysis of 
continuous linear and nonlinear systems. In this case continuous func- 
tions would be represented by sampled approximations. The analysis 
of sampled-data systems with nonlinear components offers ample pos- 
sibilities for study. These areas of sampled-data theory are not dis- 
cussed in this dissertation. 

Of the problems which are a direct outgrowth of this dissertation 
three are of particular interest: first, extention of this synthesis pro- 
cedure for systems with non-periodic operation of the various sam- 
plers in the systems and consideration of wide sampling pulses which 
cannot be replaced by impulses. Second, further extention of the 
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procedure to control the shape of the response between the sampling 
instants. The third problem concerns the need to study the zeros of 
a z-transform to develope a method of predicting the location of zeros 
from the time function. 
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CHAPTER VI 



DESCRIPTION OF THE ANALOG COMPUTER SET-UP 

An electronic analog computer study was undertaken to verify the 
synthesis methods of this dissertation. The computer was set up to 
correspond on a transfer function basis with the physical system. 

This is shown by Comparing Figure 6. 1, where a block diagram of a 
sample-data position servo is shown, and Figure 6,2 , where its com- 
puter analog is shown. The system consists of a sampler, phase lead 
network and a motor. In practice the latter device may be a commu- 
tator or any of a variety of analog to digital converters (23) (24). A 
photograph of the entire computer layout is shown in Figure 6. 3 , The 
important components of the set-up are identified on the figure. 

The computer set-up on a transfer function representation of a 
sampled data system served as more than just an aid to calculation, 
but as a direct analogy of the physical system and consequently a check 
on the analytical methods of this dissertation. Since the methods of 
this thesis are based on linear theory, the electronic analog computer 
simulation provides a more rigorous test than a physical system. No 
excuse is possible that nonlinearities obscured the results predicted 
by the theory. Also, the versatility of the computer permitted the 
testing of more possible system variations. 

As well as the system described above, a second system using a 
hold network and a third system incorporating a pure time delay as 
well as the sampler were tested. Each of the systems were designed 
to meet the same specifications and consequently the z-plane pole- zero 
configuration was the same. Two compensating networks were provided 
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Simple sample-data servo system 
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Figure 6,2. Computer set-up analyzing the system of Figure 6. 1. (all values in megohms 
and microfarads) 




figure 6.3 iLLectronic >>Jialog Computer Layout 








for each system--the first included only a phase lead network to adjust 
for rise time and peak overshoot while the second added a dipole to 
increase the velocity error constant. The details of the procedure for 
design of the compensators were explained in Section 5. 3 . Also, 
shown in Figures 5, 3 through 5.8 are copies of the output voltage ver- 
sus time curves secured from the computer set-up. These curves 
along with the tabulation of values for Table 5, 1 represent 

the observed data of the computer study. The very close check 
between the analytical results and the observed results verified the 
methods of this dissertation and no further consideration of this mat- 
ter will be made in this chapter. Rather, this chapter will limit its 
discussion to the new components which were developed to simulate 
the various system characteristics which are peculiar to a sample- 
data control system. The components referred to are (1) the Sam- 
pling unit, (2) the holding circuit, (3) the pure time delay, and (4) 
the input signal synchronizer circuit, 

(1) A front view of the Sampling unit (16) is shown on the left of 
the Figure 6, 3 . Four sets of terminals are shown since actually this 
particular device furnished four sampling switches. A top view clearly 
exhibiting the various parts making up the device is available in Figure 
6,4 . From this figure it is seen that a D. C. motor’5' drives four ganged 
differential gears through a step down gear train. The output shaft of 
each of the four differential gears drives a split cam. Adjustment of 
the cam permits control over the switch closure time. Ford automo- 
bile distributer contacts are used for the actual switch. Each differ- 



*Electric Indicator Co. , 1/25 HP, 115 volts, style ASC motor. 
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Motor Speed Control 
Potentiometer 




Fl^re 5.4 Top View of Sampling Unit 
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ential gear unit is furnished with a knob which permits displacing the 
angle alignment between its input and output shafts. Since the input 
of each differential is driven off the previous differential output (except 
for the first) this permits relative phasing between the closures of the 
four sampling switches. 

The Electronic Counter*^ shown in Figure 6,3 was used to mea- 
sure the sampling period and sampling duration. The sampling switch 
contacts were readily transferable from the computer circuit to the 
input of the counter in parallel with a battery and resistor circuit by 
means of the knife switch shown in front of the Sampling unit. This 
permitted control of the sampling period by adjustment of the potenti- 
ometer in series with the D. C, motor field. The electronic counter 
was adjusted to give time interval measurements to a tenth of a milli- 
second, A one second sampling period was chosen and it was found 
that the motor speed was sufficiently steady to hold the sampling period 

i - 

variation to less than Z. 5 milliseconds during a computer run. 

The closure time of the sampling switch was set for 59 milli- 
seconds, Employing Equation 2, 2 indicates that the actual gain of 
the computer simulated system needed to be 16, 95 times the calcu- 
lated gain in order that the pulse correspond in strength to an impulse. 
That the system would react to this pulse in approximately the same 
manner as to an impulse is assured since the principle time constant 
of the simulated system was 2, 3 seconds. The sampling period of 
one second and the 59 milliseconds pulse were chosen as a compro- 
mise between the gain and impulse approximation requirements, 

’^'Hewlett Packard 522B Electronic Counter, 
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(2) The sampler plus the holding circuit was simulated by the 
circuit depicted in Figure 6,5, It has the advantage of employing 
only one computer amplifier and the switch was furnished by the sam- 
pling unit. For proper operation the circuit functions as follows: 
when the switch closes the amplifier output voltage very quickly 
assumes the same value as the input (with reversed sign); when the 
switch opens the output voltage due to action of the capacitor will 
remain constant until the switch closes next time. It should be noted 
that the closure time of the switch should be ample to insure that the 
capacitor will charge up to its correct value of potential. With a 
switch closure time of 15 milliseconds suitable parameter values for 
R and C were found to be 0. 01 megohms and 0,015 microfa- 
rods respectively. The preceeding values give a time constant of 0. 15 
milliseconds. Since the sampling switch is closed for ten time con- 
stants it can be assumed that the difference between the input and out- 
put signal magnitudes is negligible. 

(3) Pure time delays or transportation delays can only be approx- 
imated by linear continuous networks. For this computer study the 
Fade fourth order approximant (12) was simulated by the computer 
circuit shown in Figure 6,6 (25) . The value of time delay X was 
set for 0. 5 seconds and the six operational amplifiers were secured 
from the separate Donner unit shown on the right of Figure 6. 3 . In 
Figure 6, 7 the output of this circuit to an input ramp signal is shown. 

It is significant that the output is not identically zero during the time 
delay but rather performs a low amplitude oscillation. The oscilla- 
tion is more pronounced if the input signal contains a large proportion 
of high frequencies. For this reason, the delay circuit was placed 
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Figure 6. 5. The sampler plus the holding circuit are 
simulated by this circuit 
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Figure 6. 6. Pure time delay is simulated by this computer circuit. (All values in 
megohms or microfarads. ) 
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Figure 6. 7. Input and output signals of the simulated 
time delay 
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at the output of the open- loop system permitting the system time con- 
stant to reduce the magnitude of the complementary frequency compon- 
ents introduced by the sampler. Nevertheless the output of the delay 
still contained an oscillation not present in the input signal. However, 
the system operates open circuit during this period and the signals do 
not interfere with the proper operation of the system since the delay 
assumes the proper output value by the time the sampling switch closes 
again. 

(4) Calculation of transient responses by analytical methods is 
usually performed with both the initial conditions and the input function 
zero at time zero. For a sampled data system this requires input 
functions, such as ramps, to be first applied at a sampling instant. 

For a computer study this would require that the start of a computer 
run be synchronized with a sampling instant. The circuit depicted in 
Figure 6. 8 was found to provide this synchronization. This circuit 
utilizes a second sampling switch A from the Sampling unit which must 
be synchronized with the switch used for sampling in the computer. 

The start and stop switch are both closed when a computer run is to 
be initiated. Then upon closure of the sampler contacts A , relay 
contact B closes a series circuit consisting of the coil and battery 
and thus locking the relay in the operate position while closure of relay 
contacts C starts on the computer run. The Donner computer fur- 
nishes external operate contacts which were thus employed. The com- 
puter run is terminated by placing the stop switch in the open position. 



Relay Start 




Figure 6. 8. Input signal synchronizer circuit 
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